hep-ph/0403158 



IFUP-TH/2004-2 



Probing oscillations into sterile neutrinos 
with cosmology, astrophysics and experiments 

M. Cirelli^ G. Marandella^, A. Strumia^, F. Vissani^ 

^ Physics Department, Yale University, New Haven, USA 
^ Scuola Normale Superiore and INFN, Pisa, Italia 
^ Dipartimento di Fisica dell'Universitd di Pisa, Italia 
^ INFN, Laboratori Nazionali del Gran Sasso, Italia 



Abstract 

We perform a thorough analysis of osciUation signals generated by one 
extra sterile neutrino, extending previous analyses done in simple lim- 
iting cases and including the effects of established oscillations among 
active neutrinos. We consider the following probes: solar, atmospheric, 
reactor and beam neutrinos, Big-Bang Nucleosynthesis (helium-4, deu- 
terium). Cosmic Microwave Background, Large Scale Structure, super- 
novae, neutrinos from other astrophysical sources. We find no evidence 
for a sterile neutrino in present data, identify the still allowed regions, 
and study which future experiments can best probe them: sub-MeV 
solar experiments, more precise studies of CMB or BBN, future su- 
pernova explosions, etc. We discuss how the LSND hint is strongly 
disfavoured by the constraints of (standard) cosmology. 
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1 Introduction 

High-energy colliders are the tool to discover new heavy particles with sizable couplings. New 
light particles with small couplings can be searched for in many different ways. Neutral fermions 
with eV-scale mass, called 'sterile neutrinos' in the jargon, are typically stable enough to give 
effects in cosmology, and can affect neutrino oscillation experiments.^ Both fields recently dis- 
covered new physics (in both cases around meV energies), but so far no sterile neutrinos. 

Today the most powerful cosmological probe of sterile effects is standard Big-Bang Nucleosyn- 
thesis (BBN) pp, which constrains the number of thermalized neutrinos present at T ~ 0.1 MeV 

"'^Small fermions masses are stable under quantum corrections. Cosmology, astrophysics and neutrino experi- 
ments are also sensitive to light scalars or vectors, which give different kinds of signals. 
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to be Ny = 2.5 ± 0.7. Since the uncertainty is controversial it is not clear how much Niy = 4 is 
disfavoured. 

The established solar and atmospheric neutrino anomalies seem produced by oscillations 
among the three SM neutrinos, with at most minor contributions from possible extra sterile 
neutrinos. The present 99% CL. bounds on sterile mixing, computed assuming that the initial 
active neutrino {ug} oscillates with a large Am^ into an energy-independent mixed neutrino 
cos^sl'^a) + sin^'sji^s), are sin^ 9s < 0.25 in solar oscillations (where u^. = ^e) |2| and sin^ 9s < 0.21 
in atmospheric oscillations (where z/^ = z/^) jS]- 

We relax these simplifying assumptions and study the more general 4- neutrino context. This 
demands a remarkable effort, but it is an important task for the present, for two reasons. First, 
because recent discoveries in cosmology and neutrino physics are stimulating new experiments 
that will study the new phenomena with the redundancy necessary to test the minimal models 
suggested by present data. Second, because oscillations in extra light particles are the natural 
extension of the emergent massive neutrinos scenario. 

Although we list hints that might be interpreted as sterile neutrino effects (LSND certain 
pieces of solar or atmospheric or BBN data, ...), we do not focus on any of them in particular. 
Rather, we compute experimental capabilities and constraints on sterile oscillations and compare 
them with capabilities and constraints from various present and future cosmological and astro- 
physical probes. Some important constraints are based on untested assumptions and plagued 
by systematic uncertainties. In such cases, rather than performing global fits, we identify and 
compute the key observables, trying to explain the basic physics in simple terms emphasizing 
the controversial issues, so that the reader can judge. 

The paper is organized as follows. In section |21 we briefly review theoretical motivations: as 
in the case of the known light particles, some fundamental reason would presumably be behind 
the lightness of an extra sterile neutrino. We also describe the non-standard parametrization 
of active/sterile mixing that we choose (because more convenient and intuitive than standard 
parametrizations) and describe the qualitatively different kinds of spectra on which we will focus. 

In section IHl we study sterile effects in cosmology, comparing the relative sensitivities of 
two BBN probes (the helium-4 and deuterium abundances), of Cosmic Microwave Background 
(CMB) and of Large Scale Structures (LSS). In section |3] we study sterile oscillations in solar 
(and KamLAND) neutrinos. In section El we study sterile oscillations in SN1987A and future 
supernovae. We also briefly discuss other less promising probes (relic SN background, high- 
energy cosmic neutrinos,...). In section El we study sterile oscillations in atmospheric and 
reactor neutrinos, and in short and long-baseline neutrino beams (including LSND). 

In order to avoid a unreadably long paper each section is written in a concise way and 
contains a 'Results' subsection, which can be read skipping the other more technical parts. 
From a computational point of view, exploring 4z/ oscillations is 3 -^ 4 orders of magnitude 
more demanding than usual 2 or 3u fits and therefore requires significant improvements of usual 
techniques. Subsections entitled 'technical details' describe how this was achieved. In particular 
we describe how we employ the density matrix formalism in neutrino oscillation computations, 
and explain why even in simple Su situations it gives significant advantages with respect to more 
direct equivalent approaches. 

In section [7\ we conclude by summarizing and comparing the different probes of sterile neu- 
trinos. 



2 



2 Active/sterile neutrino mixing 



We first introduce a non-standard useful parameterization of the most generic Au spectrum, and 
later present a brief review of theoretical models. 



2.1 Parameterization 

A generic 4x4 Majorana neutrino mass matrix is described by 4 masses, 6 mixing angles and 
6 CP-violating phases; 3 of them affect oscillations. 

In absence of sterile neutrinos, we denote by U the usual 3x3 mixing matrix that relates 
neutrino flavour eigenstates i^e,fi,T to active neutrino mass eigenstates 1^123 as = Uuuf {i = 
{1,2,3}, i = {e,/i,r}). The extra sterile neutrino can mix with one arbitrary combination of 
active neutrinos. 

Mixing of the sterile neutrino can be therefore fully described by a complex unit 3-versor ft 
(containing two CP- violating phases) and by one mixing angle ^s- With this parameterization 
the 4 neutrino mass eigenstates are 

{z/4 = i^s cos 9s + n£Ui sin^s q \ 

Ui = Ua[5u' - n}ne>{l - cos 9s)]iye' - sin 9sn}U^ii's 
i.e. the 4x4 neutrino mixing matrix V as that relates flavour to mass eigenstates as i^e,fi,T,s = 
V ■ 1^1,2,3,4 is 

^^n-il-cos9s)f0n sin^r?*X /f/ OX 
V — sm 9sn cos 9s J \ I J 

or, more explicitly, 



V = " ^i^iC^ " ^s) sin 9s 

Vs \ —ni sin 9s cos 9s 

In summary, n identifies which combination of active neutrinos mixes with z/g with mixing angle 

In order to understand how neutrinos oscillate in the generic case, it is convenient to focus 
on two different kinds of limiting cases, pictorially exemplified in fig. ^ 

• Mixing with a flavour eigenstate (fig. ^): n ■ u = U£ (£ = e or /i or r). The sterile 
neutrino oscillates into a well defined flavour at 3 different Am^ (which cannot all be 
smaller than the observed splittings Am'^^^^ ^^^) . 

• Mixing with a mass eigenstate (fig. ^): n ■ v = Vi {i = 1 or 2 or 3).^ The sterile 
neutrino oscillates into a neutrino of mixed flavour at a single Am^, which can be arbitrarily 
small. 

^From here on, we omit the superscript 'a' (that stands for active): it should be clear that whenever we speak 
of v^jvs mixing, this is just a short-hand for /vs- 
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Figure 1: Basic kinds of four neutrino mass spectra. Left: sterile mixing with a flavour 
eigenstate ("z/^ in the picture). Right: sterile mixing with a mass eigenstate (u2 in the picture). 

We think that our parametrization of sterile mixing, in eq. makes physics more transparent 
than other frequently used choices^. 

The oscillation probabilities among active neutrinos in the limit where the active/sterile mass 
splitting dominates, and active/active mass splittings can be neglected, are 

P[u, ^ u,) - P[u, ^ u,) - j 4|^^2||^^2j sin2(Am?4L/4E,) for £ ^ f 

and in our parametrization = sin ^s- 

Older papers studied active/sterile mixing in 2 neutrino approximation. In such a case 
6's = 7r/2 gives no oscillation effect. On the contrary, in the full 4 neutrino case 6s = 7r/2 swaps 
the sterile neutrino with one active neutrino, (e.g. t'^ in fig. ^ or U2 in fig. ^p, if 6's were there 
increased up to vr/2) affecting solar and atmospheric oscillations in an obvious way. Therefore 
large active/sterile mixing is excluded by experiments for all values of Am^4 = ml — ml (with 
one exception: the sterile neutrino mixes with a mass eigenstate z/j and the two states form a 
quasi-degenerate pair. This structure arises naturally in certain models 

In order to explore a more interesting slice of parameter space when considering sterile mixing 
with a mass eigenstate i^j, for 6^ > /t/4 we modify the spectrum of neutrino masses and replace 
{mj,m1) with {2mf — m^^mj). In such a way, the mostly active state always keeps the same 
squared mass (that we fix to its experimental value), so that in the limit = 7r/2 the sterile 
neutrino gives no effect rather than giving an already excluded effect. Physically, in our z/s/i^j 

■^Whcn studying sterile mixing with a flavour eigenstate our expression is directly related to the 'standard' 
parameterization 

V = i?34i?24-Rl4 • U23U13U12 

where Rij represents a rotation in the ij plane by angle Oij and Uij a complex rotation in the ij plane. ^14 or 
Uei gives rise to Ve/i^s mixing, 624 or [7^4 to v^/vs mixing, and ^34 or Uri to Vt jv^ mixing. 

The above 'standard' parameterization becomes inconvenient when studying mixing with a mass eigenstate. 
In such a case our parameterization is directly related to the alternative 'standard' parameterization appropriate 
for this case, 

y = U23U13U12 ■ R34R24R14 

Now 6i4 gives rise to Vi/vs mixing. Our parameterization instead is convenient because it remains simple in both 
cases. 
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plots the mostly sterile neutrino is heavier (lighter) than the mass eigenstate Vi to which it mixes 
when < 7r/4 (^s > 7r/4). When studying mixing with a flavour eigenstate we do not modify 
the spectra at 6's > vr/4 in order to obtain some other experimentally allowed configuration. For 
this reason, we restrict such plots to < 7r/4. 

We do not consider '2 + 2' neutrino mixing, namely two neutrino couples separated by a 
mass splitting much larger than ^m1^^^^^. In fact this spectrum does not reduce to active-only 
oscillations in any limiting case so that sterile effects are always sizable, and present experiments 
already exclude this possibility [Hj. When the separation among the two couples is comparable 
to Am^yj^ jjjj^, '2 + 2' is no longer a special case qualitatively different from '3 + 1'. 

We assume that active neutrinos have normal hierarchy, Am^a > 0. Finally, we assume 
^13 = 0. We verified that using ^13 ~ 0.2, the maximal value allowed by present experiments, 
leads to minor (in some cases) or no (in other cases) modifications, that we do not discuss. 
Measuring and discovering sterile effects will likely be two independent issues (however it is 
curious to note that both could first manifest as disappearance of reactor Vg)- 

These assumptions are made because we consider three active neutrinos with normal hierar- 
chy as the most plausible spectrum, and view inverted hierarchy, large and sterile neutrinos 
as possible surprises: we here study the latter one. 



2.2 Theory 

The relevant terms in the SU(2)2,-invariant effective Lagrangian that describes active neutrinos 
V together with extra light singlet fermions are 

^(^^) + ^-iyR + ^jyRLH + h.c. (3) 

H is the higgs doublet with vacuum expectation value (0,w). The first dimension-5 operator 
gives Majorana i> masses and is naturally small if lepton number is broken at a high-energy 
scale. The second term gives Majorana uji masses mjifi, and the third term Dirac ul'^r masses 
rriLR- one needs to understand why uilr and are small. 

A few theoretically favoured patterns emerge from rather general naturalness considerations. 
We consider the most generic mass matrix with LL, RR and LR mass terms. If rriLL dominates 
one obtains light sterile neutrinos with mass mg <^ ma and active/sterile mixings 6'^ ~ mg/ma. 
If mnji dominates sterile neutrinos are heavy with 6"^ ~ mg^/m^. If m^/j dominates one obtains 
quasi-Dirac neutrinos that split into couples. These 'more likely' regions can be represented as 
lines in the logarithmic (tan^ ^s, Am^^) plane, that we will use to present our results. 

A new light particle would probably be a discovery of fundamental importance, because it 
lightness is likely related to some fundamental principle, as it is the case for the known light 
particles, the photon, the neutrinos and the graviton. Attempts of guessing physics beyond the 
SM from first principles motivate a number of fermions which might have TeV^/Mpi masses and 
behave as sterile neutrinos. A few candidates are, in alphabetic order, axino, branino, dilatino, 
familino, Goldstino, Majorino, modulino, radino. These ambitious approaches so far do not give 
useful predictions on the flavour parameters in the effective Lagrangian of eq. Q. Therefore one 
needs to consider more specific ad- hoc models jT]. 

Unification of matter fermions into SO(IO) 16 or Eq 27 representations predicts extra sin- 
glets, which however generically receive GUT-scale masses. It is easy to invent ad-hoc discrete 
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or continuous symmetries that keep a fermion light [7j. One might prefer to use only ingredients 
already present in the SM. For example, the extra fermions can be forced to be light assuming 
that they are chiral under some extra gauge symmetry (that could possibly become non pertur- 
bative at some QCD-like scale, and give composite sterile neutrinos) [7]. Alternatively, the extra 
fermions may be light for the same reason why neutrinos are light in the SM j7] . Following this 
point of view up to its extreme, one can add to the SM a set of 'mirror particles', obtaining 3 ster- 
ile neutrinos [7j. In presence of Planck- suppressed corrections that mix the two sectors, mirror 
neutrinos with the same mass as SM neutrinos give rise to quasi-maximal mixing angles between 
SM neutrino mass eigenstates and sterile neutrinos, splitted by Am^ ~ rriiv'^/Mpi ~ 10~^eV 
(this kind of 1^2/1^8 effects are disfavoured by solar data). Mirror neutrinos with different masses 
from SM neutrinos can give detectable ui/us oscillation effects. The two sectors might instead 
communicate because coupled to the same heavy see-saw neutrinos: this gives massless states 
which can be mostly sterile. 

3 Sterile effects in cosmology 

Present cosmological data seem compatible with the following minimal assumption (see e.g. jH]): 
primordial perturbations are generated by minimal inflation (fiat space, Gaussian perturbations 
with fiat spectral index Ug = 1) and evolve as dictated by general relativity (with a small 
cosmological constant) and by the Standard Model (adding some unknown Cold Dark Matter). 
Global analyses performed under this assumption give precise determinations of the cosmological 
parameters and bounds on non standard properties of neutrinos. In particular, the sensitivity 
to neutrino masses and to oscillations into extra sterile neutrinos is competitive with direct 
experimental bounds. 

Present bounds are unsafe because based on assumptions which have been only partially 
tested. At the level of precision needed to derive bounds on sterile neutrinos stronger than ex- 
perimental bounds these assumptions stand practically untested. Since these bounds come from 
a few key measurements (rather than from a redundant set of different observables) compensa- 
tions among different kinds of new physics are not unnatural. For instance, the dominant BBN 
probe PP, the helium-4 abundancy jU], is plagued by controversial systematic uncertainties, and 
presently it is not clearly incompatible with a fourth thermalized sterile neutrino [TU] 

Therefore we do not try to attribute a precise probabilistic meaning to cosmological bounds 
by performing global fits. We prefer to identify and compute the key observables. This allows 
to present the underlying physics in a simple and critical way (while global fits would indirectly 
involve other cosmological parameters and data). The implications of present or future measure- 
ments can be read from our plots. Cosmology is interesting not only because gives indications 
today, but also because it will allow powerful future searches. 

BBN probes the total energy density at T ~ (0.1 ^ 1) MeV (dominantly stored in electrons, 
photons and neutrinos according to the SM) and is also directly sensitive to reactions involving 
neutrinos (e.g. UeP ^ en). Given a few input parameters (the effective number N^, of ther- 
malized relativistic species, the baryon asymmetry ns/n^ = tj, and possibly the z/^/z/£ lepton 
asymmetries) BBN successfully predicts the abundances of several light nuclei P] . Its non trivial 
success strongly indicates that primordial BBN really happened. Today 77 is best determined 
within minimal cosmology by CMB data to be 77 = (6.15 ± 0.25)10"^° jH]. Thus, neglecting the 
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lepton asymmetries (which is an excellent approximation unless they are much larger than the 
baryon asymmetry) one can use the observations of primordial abundances to test if TVj^ = 3 as 
predicted by the SM. 

Today the ^He abundancy jU] is the most sensitive probe of N^. We study also the deuterium 
abundancy, that might have brighter prospects of future improvements ^T]. For arbitrary values 
of Ny around the SM value of 3, BBN predicts P 

F„ ~ 0.248 + 0.0096 In — - + 0.013(A^'"^ - 3), (4a) 

6 15 10~^ ' \ V \ J 

^ . (2.75 ±0.13) 10- (4b) 
Fh (r//6.15 10-10)1-6 ' ^ ' 

where Yp = nA^g/nB- Both Yp and Yd are plagued by controversial systematic uncertainties. We 
do not enter into these issues and refer the reader to [121 [131 [IHl M [U] • Adopting conservative 
estimates we get 

Yp = 0.24 ± 0.01 iV^ii° ^ 2.4 ± 0.7, 

^ = (2.8 ± 0.5) 10"^ ^ ~ 3 ± 2. 
Yu 

In order to safely test and possibly rule out N^, = 4 it is necessary to either i) somewhat improve 
the determination of Yp, or to ii) improve on Id, somewhat improve on rj and on the theoretical 
uncertainty on lo- 
in our computation, we assume standard cosmology plus a sterile neutrino, and conservatively 
assume it has zero initial abundancy at T ^ MeV. Oscillations produce sterile neutrinos [T^ . 
This is a rather robust phenomenon: it is difficult to modify cosmology in order to avoid produc- 
tion of sterile neutrinos while keeping the success of BBN. In fact sterile neutrinos are dominantly 
produced at T ~ MeV, simultaneously or after neutrino decoupling. 

However, a neutrino asymmetry rj,, 8 orders of magnitude larger than t] in baryons is a 
significant extra parameter: it affects the abundancy and consequently the n/p ratio at 
freeze-out, and finally the primordial abundances. Also, a relatively large lepton asymmetry, 
around 10"^, gives extra MSW effects which can suppress active/sterile oscillations removing 
cosmological signals [T3]. More generically, allowing a non vanishing rji, any is compatible 
with the measurement of the helium-4 abundancy, just because a single measurement cannot fix 
two parameters (A'^^ and rji,). It is important to study how the situation improves when also Yd 
will be precisely measured. For example increasing N^, from 3 to 4 increases Yd by 11% and Yp 
by 5%. If the latter effect were compensated by a large neutrino asymmetry. Yd still remains 
about 8% higher. Numbers can be somewhat different, depending on how gets dynamically 
increased by sterile oscillations, which are also directly affected by a large neutrino asymmetry 
in an important way [T^ IT^ . However the general point remains. Of course one can make 
any N^, allowed by just introducing more than one free parameter; but measuring two or more 
observables would make hard to believe that new physics effects cancel among each other in all 
cases. 

In the following we stick to standard cosmology. 

For each choice of oscillation parameters, we follow the evolution with temperature of neutrino 
abundances and compute how Yp and Yd are modified. For easy of presentation, we convert their 
values into effective numbers of neutrinos, N,^^'^ and A'^^, univocally defined by the inversion of 
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eq. (01). These parameters do not necessarily lie between 3 and 4 and are employed just as a 
useful way of presenting our final results. 

We also compute another effective number of neutrinos, that parameterizes the total energy 
density in relativistic species at photon decoupling as^ 



Prelativistic P7 



l(±] 



4/3 



1 + - — N, 
8V11/ 



CMB 



(6) 



This quantity (together with other cosmological parameters) determines the pattern of fluctu- 
ations of the CMB measured by WMAP (and other experiments). Neutrinos affect CMB in 
various ways ^7]; in the future studying how neutrino free-streaming shifts the acoustic peaks 
should offer a clean way of directly counting neutrinos. Global fits at the moment imply 



jyCMB ^ 3 ± 2 (7) 

somewhat depending on which priors and on which data are included in the fit. Future data 
might start discriminating 3 from 4 neutrinos. 

Finally, neutrinos can be studied looking at distribution of galaxies because massive neutri- 
nos move without interacting, making galaxies less clustered jTH]. The effect depends on two 
parameters, that in absence of sterile neutrinos are both determined by neutrino masses: 1) the 
temperature at which neutrinos become non relativistic, ~ m^/S (so that active neutrinos 
operated when the horizon of the universe had the size that clusters of galaxies have now); 2) the 
energy density in neutrinos, Q^h^, that determines how large is the effect of neutrinos (neutrinos 
give at most a minor correction). As usual, the parameter h is iftoday/(100km/sMpc). 

Global analyses of cosmological data are usually reported as a bound on Qi,, assuming the 
standard correlation with T^, [HIE]. However, in presence of a non-thermal population of sterile 
neutrinos these two parameters are no longer universally related. E.g. there could be a little 
number of sterile neutrinos with few eV mass (heavy neutrinos affect also CMB, behaving as 
cold dark matter). This scenario has not yet been compared with present data (see [201 for 
closely related work). We assume that the bound can be approximated with the standard one 
on the following quantity:^ 

= • p] 

93.5 eV ^ ^ 

where m is the 4x4 neutrino mass matrix and p is the 4x4 neutrino density matrix, as 
discussed below. We approximate the present bound with Qi.h'^ < 0.01 (e.g. the WMAP global 
fit gives Qyh"^ < 0.76 10~^ at 95% C.L. in the standard case [H]. The bound becomes slightly 
weaker if A^^^^ = 4, or if more conservative priors or estimates of systematic uncertainties are 
adopted [IH])- The atmospheric mass splitting guarantees Qi, ^0.5 10~^ so that future attempts 
to reach this level of sensitivity are guaranteed. We assume that a 0.001 sensitivity in Q^h"^ will 
be reached. 



■^Including small effects of spectral distortions the precise SM prediction for the effective number of neutrinos 
is N^^^^ = 3.04, but we can ignore this subtlety. 

^We are neglecting CP-violation, and assuming that neutrinos and anti-neutrinos have equal density matrices. 
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3.1 Technical details 



Before discussing results, we present the main technical details. In order to compute the observ- 
ables discussed above one needs to set up the network of relevant Boltzmann equations and study 
neutrino oscillations in the early universe j2BI22]- We use our BEN code that includes all main 
effects, relying on more accurate public codes only to precisely fix the central SM values. This 
is a complicated computation because many processes proceed at T ~ MeV: active (and maybe 
sterile) neutrino oscillations, neutrino decoupling, neutron decay, electron decoupling and finally 
nucleosynthesis. The various processes have been accurately studied in the past. The extension 
from 2 neutrinos oscillations to 4 neutrinos oscillations does not involve new ingredients, but 
rewriting and implementing old ones in an appropriate way is not completely trivial due to the 
presence e.g. of oscillations at different frequencies. 

Brute force alone would not allow to explore the key observables in a vast parameter space: 
it is necessary to employ accurate approximations that one can invent understanding the physics 
of BBN. Neglecting spectral distortions possibly induced by active/sterile oscillations, the Boltz- 
mann equation for the 4x4 neutrino (and anti-neutrino) density matrix is fIT\ |22] 

P = zZh'^£ = p] - {r, (p - r)}. (9) 

where H is the Hubble constant at temperature T, z = rrie/T. The factor Z = —3 dln2;/dlns 
(where s is the entropy density) differs from 1 when the temperature of the universe does not 
decrease as the inverse of its comoving radius, namely during electron decoupling (z ~ 1) and 
when sterile neutrinos thermalize. In the flavour basis 



2E, 15 



cos^ 6'w diag {pee, Pf,i„ Ptt, 0) + 2 diag (1, 0, 0, 0) 



(10) 



We can neglect contributions from the off-diagonal elements of p j2I]. The usual MSW effect j2Sl 
gives an additional subdominant term with different sign for neutrinos and antineutrinos. The 
average over the neutrino energy spectra is performed using a Fermi-Dirac distributions, and 
therefore neglecting spectral distortions possibly caused by oscillations. 

In eq. (Q we use the standard 'anticommutator' approximation for the collision terms F, that 
describe weak ue and uu interactions that tend to thermalize neutrinos, driving their matrix 
density to its thermal equilibrium value, p*^^ = diag(l, 1, 1,0). A detailed comparison with the 
full equations [22] reveals that they are accurately mimicked by inserting the following values of 
the damping coefficients. In the equations for the off-diagonal components of p, we insert the 
total scattering rate 



Ftot ~ 3.6 forz/e and Ftot ~ 2.5 for i/^,, 

because all scatterings damp the coherent interference between different flavours. In the equa- 
tions for the diagonal components of pu, we insert the annihilation rate 



Fann-O.SG^T^ for Z/e and Fann-O.SG^T^ foiU^^r 

since annihilations are needed to change the number of neutrinos. However this procedure j221 
introduces an artificial choice of basis, giving equations which are no longer invariant under 
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rotations in the (/i, r) sector. In order to correctly maintain important coherencies among 
and Vri we introduce the distinction between Ftot and Fann in the z/^iz/,- basis (assuming maximal 
atmospheric mixing; otherwise the generalization is immediate). 

After determining neutrino evolution we can study the relative n/'p abundancy, that evolves 
according to [H [T3] 

r = zHZ^ = Tp^nil - r) - rr„ 



aZ Hn + Up 

where Tp^n is the total peue n, pe ^ nve and pve ne reaction rate (in thermal equilibrium 
the inverse process would satisfy r„^p = Fp^ne*-™'""™'''"™'''*'^^). The production of sterile neu- 
trinos affects n/p by 1) increasing the Hubble parameter H; 2) modifying the Tp^n, Tn^p 
rates, if the Ue population is depleted by oscillations. 

Finally a network of Boltzmann equations describes how electroweak, strong and electro- 
magnetic processes control the evolution of the various nuclei: p, n, D, T, '^He, ^He,. . . Rather 
than recalling here the main features of these equations, we just state (without explanation) the 
approximation we use (see also |21])- At a sufficiently low temperature T* ~ 0.08 MeV almost all 
neutrons wind up in ^He, so that its mass abundancy is given by 1^ ~ 2r(T*) with T* obtained 
solving 

180H = F 

The precise numerical value is fixed in such a way that in the SM case our simplified code 
precisely agrees with state of the art codes (that include thermal, radiative and other corrections 
corrections, smaller than the present experimental uncertainty).^ The deuterium abundancy is 
obtained with a similar technique. 



3.2 Results 

We plot the effective numbers N^, of neutrinos defined in terms of the physical observables (the 
^He and D abundances and the energy density at recombination) from eq.s Q and ©• We also 
plot the value of the present energy density in neutrinos fi^, probed by observations of Large 
Scale Structure together with CMB constraints. 

The plots have the following meaning: shaded regions have A^^"^ > 3.8 or Q^h^ > 10'^ 
and are therefore 'disfavoured' or 'excluded' (depending on how conservatively one estimates 
systematic uncertainties) within minimal cosmology. The other lines indicate the sensitivity 
that future experiments might reach. More precisely we plot contour-lines corresponding to 
K = 3.2 and 3.8 and to fi^/i^ = 10-^ and lO'^. 

recent paper |^ studied, in a four-neutrino context, how active/sterile oscillations affect the ''He abun- 
dancy. The authors of take into account spectral distortions of I's (that we neglect), and neglects other 
minor corrections (that we include), such as those related to the electron mass. They compute an 'effective BBN 
neutrino number' N^^^ , which, unlike our Nj/^° and , is not directly related to the observable helium-4 
and deuterium abundances. However, in the parameter range covered by their plots, their N^^^ should be an 
approximation to . The plots that can be compared show a reasonable level of agreement. 

At Am? ~ 10~* eV^ spectral distortions of electron neutrinos are not negligible: in the two-neutrino limit they 
make the helium-4 abundancy more sensitive to sterile oscillations (see the papers by Kirilova et al. in 
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Figure 2: Isocurves to the effective number of neutrinos produced by 2 neutrino oscillations in 
the cases Ve/v& (l^ft plot) and v^^r/i^s- Solar and atmospheric oscillations are included in the 3 
neutrino plots of fig. where the meaning of the various isolines is precisely explained. 

It is useful to start discussing the unrealistic but simple cases considered in old papers [H]. 
In fig. 121 we show the effects produced by 2 neutrino mixing: Vs/^e in fig- and z/g/i/^ or fs/^r 
mixing in fig. IJb-'^ The red dashed line shows the total number of neutrinos, A'^^^^: it essentially 
does not depend on which flavour (z/g, or Vr) mixes with z/g and is not affected by oscillations 
with Am^ < 10~^ eV^ that are too slow and start only after neutrino decoupling, when the total 
number of neutrinos is frozen.^ At this stage neutrinos can still change flavour. The difference 
between fig.|2^ and|2l3 is due to the fact that only electron neutrinos are involved in the reactions 
that control the n/p ratio. Therefore z/g/i^g oscillations that occur after neutrino freeze-out and 
that do not affect the total number of neutrinos (z/g are created by depleting Ve) affect n/p and 
consequently the ^He abundancy^ (continuous line) [Hj and, to a lesser extent, the D abundancy. 
This happens down to Am^ ~ 10^^ eV^: oscillations with Am^ < 10~^ eV^ occur after decoupling 
of electroweak scatterings, when the relative n/p abundancy is only affected by neutron decay. 

Effects are larger at 6' > 7r/4 (i.e. tan^ > 1) because this corresponds to having a mostly 
sterile state lighter than the mostly active state, giving rise to MSW resonances in neutrinos 
and anti-neutrinos (like in cosmology, also supernova z/g feel a MSW resonance for ^s>7r/4. 

^Previous papers studied the ^He abundancy and we agree with their results. We however use as a variable 
tSiV? 9 rather than sin^ 26*, so that we unify in a unique plot the non-resonant [Q < 6 < 7r/4) and the resonant 
(7r/4 <9 < 7r/2) case. 

^To be precise we should say 'the total entropy in neutrinos per comoving volume remains constant'. For 
simplicity we will adopt such loose abbreviations. 

^These region are strongly disfavoured because have a hclium-4 abundancy corresponding to Nj^'^ > 4, up to 
about 5. 
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Figure 3: Cosmological effects of sterile neutrino oscillations. We compare four different 
signals. The continuous red line refers to the "^He ahundancy (we shaded as 'disfavoured' regions 
where its value corresponds to N^, > 3.8), the violet dotted line to the deuterium ahundancy, and 
the dashed blue line to the effective number of neutrinos at recombination. We plotted isolines 
of these three signals corresponding to an effective number of neutrinos Ny = 3.2 and 3.8. The 
precise meaning of the parameter Ny in the three cases is explained in the text. The upper 
(lower) dot-dashed orange lines corresponds to Vt^h'^ = 10~^ (10~^), where fl^ is the present 
energy density in neutrinos. 
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On the contrary solar neutrinos feel a resonance for 0s<7r/4). In the past years it has been 
debated about if a neutrino asymmetry and/or large inhomogeneity develop as a consequence 
of non-linear effects, and this issue has not yet been fully clarified. Our Boltzmann equations 
assume that both these effects can be neglected. At tan ^ > 1 the bound from Q^, holds even for 
very small mixing, 6' ~ 7r/2 just because these region correspond to heavy active neutrinos. 

We now discuss how the above picture changes taking into account oscillations among active 
neutrinos. Our results are shown in fig. El the upper row refers to sterile mixing with mass 
eigenstates 1^1,2,3 and the lower row to mixing with flavour eigenstates z/g and i^^^^. 

An inspection of the upper row shows that their main features can be understood in terms 
of the (unrealistic) results in the case of 2 neutrino mixing, fig. |2l Having assumed 613 = 0, 
4z/ sterile mixing with 1/3 gives no new effects with respect to 2z/ sterile mixing with z/^^^. Due 
to solar and atmospheric oscillations, depletion due to oscillations into sterile neutrinos now 
happens in all other cases and becomes milder, because no longer confined to z/g but shared 
among all active neutrinos. Fig. shows the effects of i's/'^2 mixing (this kind of neutrino 
spectrum is plotted in fig. ^): since z/2 contains some z/g component, electron neutrinos are in 
part directly affected. Fig.|2ti' shows the effects of I's/'^i mixing: z/g depletion effects are largest in 
this last case because ui is the neutrino eigenstate with the largest z/g component. In summary: 
depletion gets transferred to all u flavours and diluted. 

Mixing with flavour eigenstates is qualitatively different, for the general reasons explained in 
section |21 We can see the effects of the solar (atmospheric) mass splitting as bumps in fig. Eli 
(jHh) where cosmological effects of z/g/z/g (z/g/z^/^.r) mixing are computed. In the case of z/g/z/g 
mixing N^, N^^^ = 4 is reached only if the sterile neutrino has a large enough Am^: the solar 
Am^y^ ^ 0.7 10~^eV^ alone is not sufficient, as also indicated by the 2u limit plotted in fig. 
On the contrary, fig. |2b shows that in the case of Ug/ t mixing a Am^ ~ Am^^^^ a; 2 10~^ eV^ 
is large enough to reach k, 4 for any value of the sterile mass. We have verified that setting 
6*13 ~ 0.2 fig.slHldo not vary in a significant way. 

3.3 Hints and anomalies: cosmology 

To conclude we list cosmological data that do not fit well into the scheme presently considered 
as standard, and that can be interpreted as manifestations of sterile neutrino effects: 

• Various determinations of the primordial helium-4 abundancy Yp try to reduce uncer- 
tainties by appropriately choosing and modeling the astrophysical systems used for the 
observation p. Some analyses find lower values of Yp corresponding to less than 3 neu- 
trinos (an effect which could come from active-sterile oscillations), even after a detailed 
examination of the systematic uncertainties. Other determinations give higher values of 
Yp, more compatible with Ny = 3. 

• X-ray cluster data seem to prefer a lower value of the parameter as than the other sets of 
CMB and LSS data, considered in [21111]. A degenerate spectrum of active neutrinos can 
alleviate this discrepancy jSHl- Otherwise, one can invoke a contribution to Q^, from sterile 
neutrinos. 
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• Decays of sterile neutrinos with mass ~ 200 MeV provide an interpretation of the reioniza- 
tion at large redshift observed by WMAP, alternative to the standard one (early formation 
of massive stars) [2ZI- 

• Decays of sterile neutrinos with mass ~ 10 MeV and abundancy = 10~^"^^ can be a non 
standard source of recently observed galactic positrons j2Hl- 



4 Sterile effects in solar (and KamLAND) neutrinos 

We compare sterile effects with present and future solar neutrino experiments and with Kam- 
LAND reactor anti-neutrino data. While computing sterile effects in reactor Uf, is straightforward, 
solar neutrinos are detected after a long trip from the center of the sun during which they can 
experience sterile effects in several different ways. This is what makes solar neutrinos a power- 
ful probe of sterile effects. In section 14.11 we describe how we precisely compute these effects. 
Section 14.21 describes how we fit data, dealing with the complication that these data contain 
a positive evidence for active/active oscillations. Results and a qualitative understanding are 
presented in section I^Hl In section we compare our results with previous analyses performed 
in limiting cases. 



4.1 Technical details 

In order to understand what happens when a sterile neutrino is added and to write a sufficiently 
fast numerical code, we must develop an analytical approximation. It is convenient to employ 
some more formalism than in the 2x2 case and study the evolution of the 4x4 neutrino 
density matrix pm written in the mass basis of instantaneous mass eigenstates. A Ue produced 
with energy E^, at radius r = vq inside the sun is described by pm = ■ diag (1, 0, 0, 0) ■ Vm 
where Vm depends on E^, and Tq. Mixing matrices in matter (Kn) and vacuum {V) are computed 
diagonalizing the Hamiltonian 



= + v^Grdiag A^^ 

2E^ 2 2 2 ^ 

and ordering eigenstates according to their eigenvalues Hi = ml^./2E^: Umi {I'mi) is the lightest 
(heaviest) neutrino mass eigenstate in matter. The evolution up to the detection point is de- 
scribed by a 4 X 4 unitary evolution matrix ^ so that at detection point the density matrix in the 
basis of flavour eigenstates i?, p = {V ■ ■ Pm{r, Ey) ■ ■ V^). where (■ ■ ■) denotes average over 
the production point. The various observables involve additional averages over neutrino energy 
and time (and consequently over different paths in the earth and in vacuum). The oscillation 
probabilities are given by P{Pe — > t'e) = Pee, -P(t'e ^s) = Pss, etc. 

We briefly describe the main steps in the computation of , focussing on the subtle ones. 
The evolution matrix can be decomposed as 

^^^arth ' '^nfacuum ' '^sun- 

Evolution in vacuum is given by "^vacuum = diag exp{—iLml./2Ey). Combined with average 
over neutrino energy it suppresses the off-diagonal element when the phase differences among 
eigenstates i and j are large. 
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Evolution in the earth is computed in mantle/core approximation, improved by using the 
average density appropriate for each trajectory as predicted in Therefore we use 'dearth = 1 
when the earth is not crossed, dearth = P ■ diag exp(— ?Lmantie"^^^./2_E',^) ■ P"!" when only the 
mantle is crossed (for a length /^mantle); and the obvious generalization when both mantle and 
core are crossed. P = VgV^ is a non-adiabaticity factor that takes into account the sharp flavour 
variation of mass eigenstates when passing from medium A (vacuum) into medium B (the earth 
mantle) . 

Evolution in the sun is more complicated because there can be various P factors at non- 
adiabatic level crossings at radii r„ (n = 1, 2, . . . , rimax), that give 



2 2 

"^sun = Pr„--- Pr2 ' diag exp(-i / ds 17^) ■ P^ ■ diag exp(-2 / ds 



The number of level crossings nmax ranges between and a few, e.g. a Uf, produced in the side 
of the sun farther from us can experience 4 crossings. When levels i and j cross in an adiabatic 
way, Pr„ = 1. If instead level crossing is fully non adiabatic Pr„ = V^(r<r„) ■ Vm{r>rn) 
is a rotation with angle a = 90° in the (ij) plane. In general the rotation angle is given by 
tan^ a = Pc/{^ — Pc), where Pc is the level crossing probability. 

So far we only presented the well known formalism [221 EHl 1^ in a non standard way appro- 
priate for applying it in numerical computations with multiple and overlapping level crossings. 
We now need to compute Pc- This last step turns out to be non trivial. By generalizing well 
known results valid in the simpler 2x2 case [201 1^ , we find that in all the parameter range Pc 
can be accurately approximated analytically (i.e. no need of numerically solving the differential 
neutrino evolution equation idp/ds = [J^,p]). Such a simple result is possible because, in the 
relevant neutrino energy range, level crossings are non adiabatic only for 9"^ ■ Am^/10~^eV^ ^ 1 
i.e. when either 6 is small or the active/sterile Am^ is much smaller than the LMA splitting. 

To compute Pc it is convenient to consider the basis of mass eigenstates in absence of 
active/sterile mixing (i.e. 6^ —>■ ot 6s ^ 7r/2, depending on which limit is closer to the value 
of 6s under examination. Here we focus on 6s —>■ 0.). This limit allows to precisely define level- 
crossings. When z/g crosses one of the active eigenstates, (a = 1,2,3) (see fig. ^ for one 
example), the level crossing probability Pc is well approximated by 

Pc = 7 = — f^ = 7-^ ^ where sm6Z = n ■ iT" sm6s. (11) 

e^-1 ' dHa/dr ' 27r|cos2^™| " s v ; 

The above equation might seem a complicated way of rewriting the well known expression for 
Pc valid in the simpler 21/ case |3n| 13 Ij. but there is one important difference. In the 2z/ case 
one can write the result in an analogous way, which contains the mixing angle in vacuum and 
the off-diagonal elements of the Hamiltonian in vacuum. On the contrary, our 7 and ^™ must be 
computed around the resonance, where Jifaa = (or around the point where adiabaticity is 
maximally violated, in cases where there is no resonance). We emphasize that reducing the full 
4x4 Hamiltonian to the effective 2x2 Hamiltonian of the 2 states that cross and computing Jifas 
is non trivial, since sterile mixing sometimes redefines the flavour of the active neutrino involved 
in the crossing.^'' In these situations it is useful to know the physical meaning of 2J^as'- it is the 

"'^'^Neglecting this subtlety would give a qualitatively wrong result e.g. in the following situation: is mixed 
with Ve and is quasi degenerate to . 
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Figure 4: We compare the usual 'active only ' global fit of solar (shaded regions in the left plot) and 
of solar plus KamLAND data (shaded regions in the right plot), with the Gaussian approximation 
employed in this paper (dashed lines) 

minimal difference between the eigenvalues of the two states that cross. In order to elucidate 
the physical meaning of the crossing angle we emphasize that it can also be extracted from 
the scalar product between the flavour vectors of the two matter eigenstates i and i + 1 that 
cross: sin 9"^^ = v*^^ (r < r„) ■ z/^^ (r > r„) and cos 6^™ = v*^^ (r < r„) ■ Vrn,+^ (r > r„,). 

When the active/sterile mixing is large, Pc is sizable only if Am^^ <^ AmLMA- ^^is case 
our expression (fTT|) reduces to the standard 2-neutrino formula. In the narrow resonance limit 
(e.g. for 6'2^ <^ 1), eq. (jlip reduces to a Landau-Zener form Pc — e~'^^l'^: our expression for the 
adiabaticity factor 7 holds for a generic 2x2 Hamiltonian. In the sun, resonances with close 
to 7r/4 happen only in the quasi- vacuum region at the border of the sun. This will be no longer 
true when studying supernova neutrinos, that will necessitate an extension of eq. 

4.2 Experimental data and fit procedure 

We fit all available latest data: 

• The SNO CC+NC+ES spectra |22]5 divided in 34 bins (17 energy bins times 2 day/night 
bins). 

• The total CC, NC and ES rates measured by SNO with enhanced NC sensitivity 

• The Super-Kamiokande ES spectra Pl], divided in 44 zenith-angle and energy bins. 

• The Gallium rate [32], -Rca = (68.0 ±3.8) SNU, obtained averaging the most recent SAGE, 
Gallex and GNO data. 
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The Chlorine rate EE], i?ci = (2.56 ± 0.23) SNU. 



• The KamLAND reactor anti-neutrino data, divided in 13 energy bins with prompt energy 
higher than 2.6 MeV |27I. 

We revised solar model predictions and uncertainties |3S1 including the recent measurement of 
the ^^N(p, 7)^^0 nuclear cross section [22], which reduces the predicted CNO fluxes by roughly 
50%. 

Data are compared with predictions forming a that takes into account statistical, system- 
atic and theoretical uncertainties (on the total solar neutrino fluxes and on the spectrum) 
and their correlations [21] . We plot the ^is function of the 2 parameters that describe sterile 
oscillations, marginalizing the full with respect to all other sources of uncertainty including 
the LMA parameters Am"^^^ and ^sun. This step is of course not performed in usual analyses. 
Proceeding in a fully numerical way, it would be too demanding for present computers. We 
can however approximate all observables with a first order Taylor expansion around the best-fit 
LMA point, since experiments allow only relatively minor shifts from it. In this way marginal- 
ization over Am^2 ^^"^ ^sun can be performed analytically, using the same Gaussian techniques 
commonly employed for other 'systematic' parameters. Fig.a|3] shows that performing this lin- 
earization (dashed lines) we obtain a satisfactory approximation to the usual active-only global 
fit of solar neutrinos (continuous lines). With more data this approximation will become more 
and more accurate. 

Solar u and reactor u data cannot be fit independently, since both sets of data depend on 
solar oscillation parameters. Inclusion of KamLAND data presents a slight complication: due to 
poor statistics, the Poissonian distribution must be used. Nevertheless small deviations from the 
LMA best fit (due to sterile effects and to systematic uncertainties) can be taken into account 
in Gaussian approximation, obtaining 

= 2(t - r In t) ~ 2{to - r In to) + 2e(to - r) + re^ 

where r are the observed rates and t = to{l + e) the predicted rates. In our analysis e is the 
correction due to sterile neutrinos. Fig.|3jD shows that we can accurately analytically approximate 
the usual active-only global fit of solar plus KamLAND data. 

4.3 Results 

We start recalling how LMA oscillations behave in absence of sterile neutrinos [221. -^^S- Efe^ 
shows the composition of the two neutrino mass eigenstates. For sake of illustration we assumed 
normal hierarchy, mf <^ m\ ~ Aml^^ ^ m\ ^ Aml^^. At higher neutrino energies, 

E,->E,^ Aml^JGpN^'''' ~ few MeV, 

matter effects dominate around the center of the sun, r < 0.2i?sun where neutrinos are produced 
as I'e ~ t'2m. The LMA level crossing at r «i 0.2i?sun is adiabatic (in fact the solar mixing angle 
^sun is so large that fig. [Jl does not look like a level crossing) so that neutrinos produced as z/2m 
exit from the sun as 1^2 = sin^sun + cos^sun '^^i,t i.e. Pee = sin^6'sun. This limit roughly holds 
at energies probed by SNO and SK: e.g. their total rate is Pee ~ 1.15 sin^ 6'sun. At Ey<E^ Pee 
increases and vim contains some Ve because matter effects are no longer dominant. At i?^ <^ E^ 
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Figure 5: Level crossing schemes. The right plot shows qualitatively the effective anti-neutrino 
masses in a supernova. The left plot shows the effective neutrino masses in the sun. We as- 
sumed hierarchical active neutrinos (i.e. mi = and m2 = (Aml^^J^^"^ ) and plotted the two 
matter eigenstates that give rise to LMA oscillations for two different values of E^: lOMeV and 
0.86 MeV, the energy of the main Beryllium line. Colors indicate the flavour composition. An 
extra sterile neutrino with small mixing is represented by an horizontal line with height equal to 
its mass. 



matter effects are negligible and one gets averaged vacuum oscillations, P^e = 1 ~ | sin^ 2^sun- 
This energy range has been explored by Gallium experiments. 

We now discuss how to understand qualitatively the sterile/active mixing effects If 
sterile/ active mixing is small, the mostly state is represented by adding one quasi-horizontal 
line to fig. E^. Depending on its height (determined by the mass of the sterile neutrino) the 
mostly sterile level crosses one or none of the two mostly active neutrinos (for all relevant neutrino 
energies the sterile state does not cross both active neutrinos) after or before the LMA resonance 
(or after and before). In each case one can understand the behavior of the survival probabilities 
from the level-crossing scheme: in the example plotted in fig.l^ a neutrino produced at tq ^ 0.2 
experiences a single level crossing at vi 0.3. This example corresponds to the case considered 
in jin]: a sterile neutrino weakly mixed with Ui, with mass splitting Aml^ somewhat smaller 
than Am^2- ^ solar z/g produced at r ~ 0.2i?sun contains a uim. component that crosses the Us 
state once, getting partially converted into Ug. This gives a dip in the survival probability at 
intermediate energies: at low energies one has averaged vacuum oscillations (negligibly affected 
by the small sterile mixing angle), at high energies matter effects dominate so that — V2m that 
does not cross the sterile level. Fig. [7|A.,B show examples of this behavior. Even sticking to the 
case of Vs/ mixing, qualitatively different effects are present for other values of the oscillation 
parameters, studied in fig. IHK- E-g- the example in fig. [7p illustrates the case discussed in |H1[7]: 
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Figure 6: Sterile mixing: effects in solar neutrinos. No statistically significant evidence 
is found. Shaded regions: excluded at 90, 99% C.L. Coloured lines are iso-curves of a few 
promising signals. Continuous red line: A^^^ differs from LMA by 0.005. Continuous blue line: 
0.02 day/night asymmetry at Borexino. Dashed violet line: seasonal variation at Borexino with 
0.02 amplitude. Short-dashed green line: P^e at sub-MeV energies differs from LMA by 0.02. 
Other signals are discussed in the text. The letters A,. . . ,F indicate sample points, studied in 
detail in fig. 
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Figure 7: A few samples of still allowed sterile effects in solar neutrinos. We plot, as function 
of the neutrino energy, P{ve ^e) (decreasing red curve), P{ve ^ti,T) (increasing blue curve) 
and P^i^e ~* ^s) (lower black curve). The continuous (dotted) curve are the values during day 
(night). The sample points A,. . . , F are drawn in fig. as dots. 

large and small Aml^ ~ 10~^^eV^. In all these examples sterile neutrinos manifest at low 
energy, E^<E^:. There is a general reason for this behavior: in absence of sterile effects, only 
at such energies LMA oscillations allow a ui component in the solar neutrino flux. 

We also consider a sterile neutrino mixed with 1/2 or z/3. Us/ui mixing and z/s/z^2 mixing affect 
solar neutrinos in similar ways. In fact the z/]" and neutrino eigenstates in matter both 
typically contain significant ui and 1^2 components. There are some differences in 0{1) factors, 
which produce the difference between fig. IHt {i^s/^i) and b The most evident difference 

in the shape of the excluded region is due to our choice of the parameterization: when studying 
Vs/v2 mixing we produce plots with Am\^ on the vertical axis, so that small values of Am\^ 
correspond to quasi-degenerate with V2 (rather than with vi). 

A sterile neutrino mixed with 1/3 gives much smaller effects in solar neutrinos, so that we do 
not show the corresponding plot. This happens because matter effects negligibly mix 1/3 with 
Vim or V2mi SO that MSW resonances are highly non adiabatic i.e. ineffective.^^ 

As discussed in section |21 sterile mixing with one active mass eigenstate is a special configu- 
ration. Therefore we also consider a sterile neutrino mixed with a flavour eigenstate: Vei '^n ot 
Vt-. In such a case there are active/sterile oscillations at multiple Am^ values, which (in view 
of the observed mass differences among active neutrinos) cannot be all small. This is the main 
difference with respect to the previous case, and implies that sterile oscillation effects are present 
even for Am?4 = 0. To understand better this point let us consider the case Am^4 — > 0. In this 

^^Narrow dips in Pee{E^) are possible close to specific energies such that the mostly sterile state crosses active 
neutrinos when their position-dependent eigenvalues are maximal. A look at fig. |Sla, might help to understand 
this issue. 
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limit neutrinos exit from the sun as an incoherent mixture of z/2 and z/i, which both contain some 
sterile component. At E^, ^ E^, LMA is fully effective and neutrinos exit as pure U2, so that 
vacuum oscillations related to the small Am^4 have no effect. On the contrary, at Ei^ < E^, there 
is a Ui component which experiences vacuum oscillations with the mostly sterile state. These 
vacuum oscillations affect the flux (if z/g is mixed with Ug as in fig. IHt), or the u^^r flux (if Us is 
mixed with i/^ ,- as in fig. Etl). In general, a mostly sterile state significantly mixed with z/i and 
almost degenerate to it, Aml^ ~ 10~^°eV^, gives significant spectral distortions, but only below 
the energy threshold of SK and SNO. This is illustrated in flg. [7|by the sample points D, E and 
F. 

Fig. studies Vs/i^e oscillations for generic values of Am^4, and is similar to flg. int,b apart 
from the difference discussed above. Vs/i^ii mixing affects solar neutrinos in the same way as 
Vs/t-'t mixing, and both cases give relatively mild effects, as shown in fig. IHti. 

It is interesting to notice that there are sizable and specific earth matter effects if the sterile 
neutrino has a large mixing e.g. with Vf,, and is quasi-degenerate e.g. to Vi. In fact, the two 
quasi-degenerate neutrinos have different amounts of active and sterile components; due to 
matter effects, in the earth the two states are no longer quasi-degenerate, giving oscillations 
with wave-length \/27i/GpN^ ~ 6000 km, comparable to the size of the earth 

In all cases fig.slHlshow the regions already excluded at 90 and 99% C.L. (2 dof) by present 
data (these regions are precisely defined by > Xmin + 4-6 and 9.2).^^ We emphasize that 
most of the excluded regions are disfavoured only by combining high-energy (SNO, SK) with 
low-energy (Gallium) solar data: each kind of data-sets alone tests only a minor region. This 
means that present data are already able of fixing oscillation parameters with some redundancy, 
partially testing the LMA hypothesis. Only little regions are excluded if we fit only high-energy 
data, or only low-energy data. In fact, sterile oscillations that crucially involve ui affect solar 
neutrinos only at Ei, < -E"* because only in this energy range the sun emits Ui. 

The other contours in fig. sIHl show which regions can be explored by some future experiments. 

• In view of the previous comment it should be not surprising that the relatively more 
powerful future measurement (dotted green line) is an improved measurement of Pee at 
sub-MeV energies. Taking into account solar model uncertainties and plans of future 
experiments, we assumed — maybe optimistically — that it will be possible to see 0.02 
shifts from the LMA value of Pee-^^ 

• The red continuous line shows what can be achieved by measuring the day/night asymme- 
try (normalized as 2(i?day — -Rnight)/(-Rday + -Rnight)) at Boron energies with 0.01 precision 
in a future Mton water Cerenkov experiment. We have considered this test because the 
needed experiment seems highly motivated by various other considerations. 

-'^^Using instead > Xlma + ^-^ ^^'^ 9.2 would give slightly weaker constraints. One might think that this 
latter procedure should be preferred because gives more 'robust' bounds. It gives instead under-constraining 
bounds. More robust bounds should be obtained demanding a higher C.L. to the correct statistical test, rather 
than inventing 'more robust' tests. A correct 90% C.L. bound does not hold with more than 90% probability. 

^■^As in |42j we averaged Pee assuming that pp neutrinos will be detected by elastic scattering on electrons; 
similar results hold for other possible techniques. Assuming that LMA is the end of the story, ref. 021 found 
that feasible sub-MeV measurements would not have a significant impact on the determination of the oscilla- 
tion parameters. We here find that sub-MeV measurements are instead crucially important for probing sterile 
oscillations. 
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Table 1: Rough classification of possible sterile MSW resonances in the sun. rris, mi andm2 are 
respectively the masses of the mostly Ug, t^i o,nd states; E^ ~ few MeV is the LMA critical 
energy and 6'ad ~ (Am^/10~^ eV^)^/^. No effect is present if 9s -C ^ad- 

• Coming to near-future experiments, the other lines show what Borexino (and/or possibly 
KamLAND) can do by studying (mainly) Beryllium neutrinos. We do not show the impact 
of a measurement of the total rate, which has a relatively large theoretical error [211^'^ and 
instead focus on signals that LMA predicts to be unobservable. The region inside the 
continuous blue line has a day/night asymmetry in the Beryllium rate larger than 0.02. 
In the region inside the dot-dashed blue line the Beryllium rate shows anomalous seasonal 
variations with amplitude larger than 0.02. The physics behind these effects is similar 
to the one well known from discussions of LOW and (Q)VO z/g u^j^^r oscillations (now 
excluded) . 

• We do not show results for a few other signals, that seem less promising than the ones 
discussed above. One can measure better the energy spectrum of Boron neutrinos, measure 
Beryllium neutrinos both in NC and CC scatterings, and look for day/night or seasonal 
variations in pp neutrino rates. Active/sterile oscillations with Am^ ~ 10~^ and tan^ 6s ~ 
0.1 can distort the i^e energy spectrum in KamLAND or in future reactor experiments. 

In table ^ we qualitatively classify how MSW resonances with a sterile neutrino affect solar 
neutrinos. The main variables are: which state is crossed, how large is the mixing. Fig. [7| shows 
a few examples of specific oscillation patterns compatible with present data. 

4.4 Hints and anomalies: solar data 

Finally, we have studied if present solar data contain hints of sterile neutrino effects. There are 
two ways of searching for new effects: data-driven or theory-driven. 

In the data-driven approach one performs a goodness-of-fit test, that should reveal if data 
contain some generic indication for new physics beyond LMA oscillations. In practice, this means 
carefully looking if data contain hints of anomalous results. Present neutrino data contain the 
following non-statistically-significant hints: 

i^Por values of the oscillation parameters around the 'vertex' of the MSW 'triangle' (at Anil^ « 10^^'^ eV^ 
and 6s ~ 10^^) sterile neutrinos can give a deep dip in Pee at Ebc = 0.863 MeV. Such effects can significantly 
reduce the total rate at Borexino compatibly with present data and without giving different kind of signals, such 
as day/night or seasonal variations. 

^^In the past, global analyses of solar data reported the result of a Person test. Today it would tell that 
even the LOW solution (which has been excluded) is perfectly acceptable. Therefore the fact that, according to 
the test, LMA gives a good fit is not a useful information. In fact, as discussed in the test is only 
sensitive to variations in the x^ larger than v^, where N ~ 100 is the number of data points. In the present 
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• The rate measured by the Chlorine experiment is 1.4 experimental standard deviations 
lower than the best-fit value predicted by active-only oscillations.^^ 

• The lack of an up-turn in present z/g energy spectra does not give a statistically significative 
hint for new effects additional to LMA oscillations. To get a feeling of how accurately SNO 
and (mainly) SK data test the energy spectrum around ~ lOMeV we fit these data 
assuming 

P{i^e Ve) = sin' ^sun + a{MeV / = 1 - P{ue v^^r). 

At energies probed by SK and SNO LMA oscillations have this form with a ~ 2. Data 
give a = 1 ± 2.5. 

• A 2 ^ 3cr indication for sterile neutrinos might appear if analyses that suggest higher 
values of the nuclear factors 5*17(0) and 5*34(0) (and therefore a Boron flux higher than 
what measured by SNO) were confirmed 

• According to some reanalyses, solar neutrino rates of various experiments show some time 
modulation at specific frequencies. However subsequent analyses performed by some ex- 
perimental collaborations do not corroborate these claims |13]. 

The data-driven approach cannot see new physics that manifests giving small corrections to 
many observables. These diffuse minor effects can be seen if one knows what one is looking 
for. More precisely, one performs a fit assuming a theory with n extra parameters and looks if 
the best-fit improves by more than n. Our plots in fig.s IHl explore a few n = 2 representative 
slices of the full 4-dimensional parameter space of relevant active/sterile oscillation parameters: 
various regions (not shown) are favoured in a non statistically significant way. The best fit 
values, indicated by the green dots in fig.s IHl have a lower than best-fit LMA oscillations at 
most by 3 units. In order to study the general 4z/ case we also computed local minima of the 
using numerical minimization techniques. We have not found statistically significant indications 
for an extra sterile neutrino. 



4.5 What is the bound on the sterile fraction in solar oscillations? 

Assuming that solar neutrinos oscillate into v,, — > sin a z/g -|- cos a u^^^r the question is: what is the 
bound on the 'sterile fraction' rjs = sin' a? We here compare with previous results. Oscillation 
effects at Boron energies can be parameterized as 

=$BPee, =$B(l-Pee)(l-r/.), = $B ( 1 " Pee)r/s. (12) 



SK and SNO have measured and the total flux of active neutrinos ^ue,^,T- these two mea- 
surements alone cannot determine the three unknown quantities Pee, t]s and $b (unoscillated 
total Boron flux). Adding the solar model prediction for $b gives 

Vs^^ ^^0±0.2 (13) 



case N ^ 1: this explains why the test is so inefBcient and one has to perform less standard more efficient 
tests. 

^^Beforc revising solar model predictions including the recent LUNA data [21], the Chlorine rate was lower by 
1.7(7 (according to our results, in agreement with other similar analyses, see e.g. ^Ol)- The inclusion of LUNA 
results negligibly affects the global fit of solar and KamLAND data in terms of active/active oscillations. 
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where ?7s is small because $b agrees with the measured total t'e.^t.r flux. In line of principle this 
agreement might be accidental: $b could be larger than what predicted by solar models, and a 
fraction of it could oscillate into sterile neutrinos leaving a reduced ^ ^ . In view of all other 
constraints (on spectral distortions, Gallium rates,. . . ) it seems difficult to realize this scenario 
with oscillations. 

Most previous analyses prefer to obtain constraints on ?7s only from this second argument. 
They add one sterile neutrino mixed only with u^ ^- and separated by a large mass splitting: 
Am^4 ^ Am^yj^. This special configuration gives an energy- independent ?7s and an energy 
dependence of Pee similar to what obtained for ?7s = 0. This choice does not allow to fit low 
energy Gallium data together with high energy SNO and SK data: the resulting bound on ?7s 
is 

r/s^0±0.1, (14) 

somewhat stronger than the bound in eq. (fT^ obtained using only the solar model prediction for 
the Boron fiux. We stress that this constraint on rj^ can be relaxed by making milder assumptions 
on the active/sterile oscillation parameters. For example, even still maintaining Am? ^ Am"^^^ 
r)s is energy-dependent if the heavy Ug mixes with Ue, or with the mass eigenstate ui {ris is about 
2 times larger at lower energy, see fig. Ej) or with i>2 (r/s is about 2 times larger at higher energy). 

When we focus on the special configuration considered by previous analyses we get x^iVs = 
1/4) — x^iVs = 0) = 8.3 (including in the fit the BPOO prediction for the total Boron fiux) and 
6.9 (not including the prediction for the Boron fiux: this value of the Ax^ agrees with other 
analyses [2j)- In fact, in the special configuration our parameter Og is related to the sterile 
fraction as sin^ 6's = 2r]s (for maximal atmospheric mixing). We produced fig. IHl including in our 
data-set the solar model prediction for the total Boron fiux. Dropping it would give only minor 
modifications, as a comparison of eq. (fT^ with eq. (fTH) indicates. 

5 Sterile effects in supernovae and other astrophysical 
sources of neutrinos 

Having discussed sterile effects in the sun, it is useful to discuss the analogous effects in core- 
collapse supernova (SN) focussing on the differences between the two cases and emphasizing 
what one loses and what one gains studying SN neutrinos. 

• Present detectors can only study SN neutrinos from our galaxy, and possibly from nearby 
galaxies: such SN neutrinos have a duty time of C(10) seconds every C(IO^) seconds jlHl 
EH] . This makes backgrounds less problematic, but allowed to detect so far only 0(10) 
SN1987A events [52] • Running solar neutrino experiments could detect thousands of events 
from a future SN exploding at distance D ~ 10 kpc. An even more impressive harvest 
of data could come from a future Mton water-Cerenkov detector or from other more SN- 
oriented future projects |T7] . 

• In a SN the matter density grows from zero up to nuclear density (p ~ lO^'^gcm"^) in 
the stiff inner core: matter effects are important for all the mass range that we consider 
(up to about 10^ eV). In this mass range, and given the typical SN neutrino energy 
of ~ lOMeV, active/sterile MSW resonances occur outside the neutrino-sphere (roughly 
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defined as the regions after wliich neutrinos freely stream, p <^ 10^^ gcm""^). In particular, 
conversions inside the neutrino-spheres (which could drain almost all the SN energy) do 
not take place. We only receive neutrinos that come out of the neutrino-spheres on the 
side of the SN closer to us. Unlike the case of the sun, it is not necessary to average 
over the production point, because the production and the oscillation regions are spatially 
separated. 

• While only Uf. are produced in the sun, SN produce all active u and roughly in similar 
amounts [IHl- They mix and convert among themselves, and possibly with sterile neutri- 
nos 1^ 1511 15^ I52j . Present experiments can accurately study z/g. In fact in the energy 
range me <^ <^ rrip relevant for SN neutrinos, UeP — > en scatterings allows to detect Ue 
and to measure their energy. Detecting other neutrinos (e.g. via ue scattering, deuterium 
dissociation, z/g absorption on carbon) is possible and future projects could give significant 
information. In the following, we focus on the Ue flux. 

• ^e,/i,T experience no active/active matter resonance (unlike i^e,^i,T] we are assuming normal 
hierarchy, and equal initial fluxes of muon and tau- neutrinos): their effective Am^ in 
matter increases in a monotonous way when the matter density grows. The sterile effect 
that can more strongly affect the Ue rate is a low-density Us/ vi MSW resonance, possible 
when the mostly sterile neutrino is lighter than ui (i.e. > 7r/4 in our parameterization). 

Furthermore, due to the peculiar composition of the inner part of the mantle (deleptonized 
matter) the Ve potential changes its sign in the deep region of the mantle, and it does so 
in a very steep manner |^ EHl HH ED] (see fig. Eb)? adding a MSW resonance also for 
9s < vr/4. Although details are uncertain, this is a robust prediction. It implies that 
always meets a "sharp" resonance with z/e in the deep region of the mantle I52j. 

• SN neutrinos are emitted from galactic distances, allowing to probe vacuum oscillations 
with Am^ as low as 10~^® eV^ (the precise value depend on the SN and on its distance 
from the earth: different SN probe different ranges). 

SN neutrino reach energies higher than solar perhaps up to 100 MeV, and could therefore 
probe earth matter effects and spectral distortions. 

• While the sun is essentially static, an exploding SN is a dynamical environment (neutrino 
light-curve evolution, passage of the shock wave...); including the time dependence in the 
neutrino fluxes and in the matter density profile is too much demanding and probably 
useless for our purposes, given the poor knowledge of the details. We focus on a typical 
SN configuration, which includes all the characteristic features of the SN cooling phase. 

A reliable prediction of the emitted fluxes is still lacking and the explosion mechanism is not 
yet under control. However, it is relatively easier to predict the fe.^i.r, ^e,^^,T energy spectra, 
which in thermal approximation only depend on their cross sections. On the contrary their 
total fluxes also depend on the SN density profile; in particular the thermal approximation 
does not imply equipartition among different flavours. The underlying complexity of SN 

"'^^The resonances with the sterile state would enter in the neutrino-spheres (in the inner core) for Am^ > 10^ eV^ 
( > 10'' eV^ respectively). 
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explosions makes hard to reduce theoretical uncertainties, posing a threat on the usefulness 
of precise SN v experiments as tools for studying oscillations. 

As a last warning, we have to recall the reluctance of the neutrino data from the only SN 
we know to fit into a simple and straightforward interpretation. The main reasons are: 
the energy distribution of Kamiokande and 1MB look different; both experiments have an 
excess of forward events; the 5 LSD events cannot be accounted for. 

As for the future, there are various interesting observables, and it is difficult to guess on which 
signals we should focus. This will probably depend on compromises among future experimental 
and theoretical capabilities. 

Several SN- related bounds on the active/sterile neutrino mixing have been considered in the 
literature (SHESlEllEni- Focussing on the cooling phase neutrinos, they include (i) the detection 
of Vf. from SN1987a in the Kamiokande and 1MB experiments [Snj, which sets a constraint on 
the portion of neutrinos that oscillate into (ii) r-process nucleosynthesis: a fraction of the 
heavy elements in nature is supposed to be synthesized in the region surrounding the core of 
the exploding stars, provided that the electron fraction Y^. < 0.5; even small modifications of 
neutrino fiuxes affect the process, so that the request of a successful nucleosynthesis has been 
used to set limits (see e.g. [S31ini]); (iii) re-heating of the shock: in the delayed-shock/neutrino- 
driven picture of SN explosion, the flux of neutrinos and antineutrinos from the early stages of 
the accretion phase (which should carry 10 -=-20 % of the total flux) are responsible for the actual 
explosion of the star, pushing from below the stalling shock wave; since and z/g are the most 
effective in this (interacting with charged and neutral currents with baryonic matter), a cut of 
their flux would prevent this mechanism from working [^EHEZl; (iv)... We will only consider 
the bound from direct observation, which we consider robust enough.^® 

5.1 Technical details 

The analysis of SN neutrinos can be carried on in a way similar to that of solar neutrinos, after 
adapting the computational procedure to SN peculiarities, which introduce some complications 
and allow some simpliflcations. We must follow the fate of the neutrinos emitted from neutrino- 
spheres along their travel through the star matter, the vacuum and the earth. The density matrix 
formalism, already described in the 'solar' section 14. ![ automatically handles the complication 
that the SN initial neutrino flux has a mixed flavour composition. In the SN case one should 
follow the evolution of two 4x4 density matrices: p{E^) for neutrinos, and p{Ey) for anti- 
neutrinos. We focus on the total Vf, rate, as measured by UeV scatterings with the cuts and 
efficiency of the Kamiokandell experiment. The cross section is taken from |58j . 

MSW resonances affect the neutrino density matrix, possibly introducing off-diagonal ele- 
ments which are however averaged to zero by large oscillation phases. Therefore we can combine 
probabilities rather than quantum amplitudes: the 4x4 density matrix can be replaced by 

^^Indced, for instance, alternative sites for effective nucleosynthesis have been repeatedly proposed. Moreover, 
the position of such a bound could depend quite heavily on the complicated details of SN dynamics and, in 
addition, on the interplay of sterile oscillations with it, see below. On the same footing, since the details and the 
nature itself of the neutrino-driven explosion mechanism are still to be fully understood, we do not consider it 
here as a robust constraint, not mentioning that it would require a demanding understanding and simulation of 
the evolution of the SN mantle. 
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a 4-vector containing its diagonal elements, $ = (pn, P22, P33, P44) and the evolution equation 
p = ^ • p° • = . . . by 



$ = ^„---^i-$° where {^n)ij = \{'^n)ij? (15) 

are 4x4 matrices of conversion probabilities. At the production region, matter effects are dom- 
inant, so that matter eigenstates coincide with flavour eigenstates (up to a trivial permutation): 
$° = ($p_^, $p^, $p^). $ are the fluxes of neutrino vacuum eigenstates reaching the earth 
surface, related to the fluxes of flavour eigenstates by an incoherent sum weighted by the the 
neutrino mixing matrix V: 

($p,,$p,,$p,,$p,)^=5:iK,p$, (16) 



If neutrinos cross the earth, earth matter effects [59] may reintroduce coherencies among the 
various fluxes, and one must return to the general quantum expression. Similarly, if two states 
have Am^ < 10~^^ eV^ vacuum oscillations do not give large phases: evolution in the outer region 
of the SN and in vacuum must be described keeping the off-diagonal components of the matrix 
density. In practice, we found convenient to evolve the density matrix in the instantaneous mass 
eigenstate basis, setting to zero its off diagonal elements when they accumulate large phases. 
This procedure is described in greater detail in the 'atmospheric' section lUTTl 

The effect of active/active oscillations is well known the LMA resonance is adiabatic 
and partially swaps active neutrinos, giving^^ 



= COS^ ^sun$™ + Sin^ ^sun*™ = ^sun$™ + COS^ ^sun*™ (17) 

Active/sterile MSW resonances can occur after of before the LMA resonance. 

Concerning the initial fluxes, the accurate results of simulations are usually empirically ap- 
proximated by a so called "pinched" Fermi-Dirac spectrum for each flavor a = z/g, z/e, t^x {^x 
collectively denotes f/x,r, i^^i.t) inOlllE] 

■{E^,t} = N{ria) 



dtdE^^ ' "°'77r4T4e^-/^-'?- + 1 ' ' 

where the pinching parameter 77^ takes the typical values ?7i/<, ~ 5 — 3, 77^^ ~ 2.5 — 2, ?7i^^ ~ — 2. 
For ?7a = the normalization factor is = 1 and T = (i?j,)/3.15. Based on the recent results 
of jOI], we adopt the following average energies and total luminosities for the various neutrino 
components at the time of the snapshot of fig. |H1 (see below) 

(E,J~12MeV, (EpJ~14MeV, ~ 14MeV 

L^^ ~ 30 ■ 10^1 erg sec'S Lp^ ~ 30 ■ 10^^ erg sec'S L^^ ~ 20 ■ 10^^ erg sec'^. ^ ' 

In accordance with numerical calculations, we shall assume that the ratios of luminosities do 
not vary much during the whole emission. The initial flux of sterile neutrinos is assumed to 
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Figure 8: Supernova profile. Fig. Density p{r) and electron fraction Yf,{r) from \5(J^ . 
Fig. O).- matter potentials in the SN, for i/^ (red solid line) and for u^^r (blue solid line). The 
dashed lines are the analytic modelization that we adopt. 



be vanishing, as a consequence of tlie fact that matter oscillations only take place out of the 
neutrinosphere. 

Neutrinos traveling in the SN matter experience the MSW potentials jlHl 1^ 

K = v^GpnB (31; - 1) /2, Vr = V^ + V^r, . . 

= v^GpriB (i; - 1) A V, = 0, ^ ^ 

where is the baryon number density {ub = p/friN where ^ 939 Me V is the nucleon 
mass^^) and Y^. = {N^- — Nf,+ )/nB is the electron fraction per baryon. Antineutrinos experience 
the same potentials with opposite sign. 

The difference V"^^ in the i/^ and Vj. potentials, which appears at one loop level due to the 

^^We are assuming normal hierarchy and Oij, = 0. The situation becomes more comphcated if instead 9i3 > 1°, 
because the atmospheric resonance starts changing the resuh in eq. H17|) . For 1° ^ ^13 ^ 1 it is adiabatic. If 
Am|3 > (normal hierarchy) it affects only neutrinos by giving a full conversion, ^ff^°^ = If instead 

Am23 > (inverted hierarchy) it affects only anti-neutrinos, by giving a full conversion. <f>|^*'"^ = 

^°To be precise, shortly after the collapse the electron neutrino matter potential in the very center of the 
core is positive, due to the contribution of the trapped neutrinos themselves. This configuration lasts for a short 
transient period, until the neutrino diffusion depletes their abundance and carries the potential to negative values, 
where it stays. A small fraction of the electron (anti) neutrinos produced in the deep core could then oscillate 
into sterile states and constitute a non vanishing flux injected in the mantle. 

■^^For the very high densities closer to the core, ttln should be replaced by a (quite different) effective nucleon 
mass. We can neglect this refinement in the regions of our interest. 
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different masses of the muon and tau leptons is, according to the SM 




(21) 



The effect is not irrelevant in the inner dense regions: for densities above p ~ 10^ gcni"'^, the /ir 
vacuum mixing is suppressed. 

A crucial point concerns the characteristic of the matter density and of the electron fraction in 
the mantle of the star. We adopt the profiles represented in fig. |H1 and we model them with 
analytic functions that preserve their main features. Namely, the density profile decreases 
according to a power law out of the ~ 10 km inner core (where instead it has a roughly 
constant, nuclear density value). At much larger distances the density profile gets modified in 
a time-dependent way by the passage of the shock wave. Present simulations have difficulties 
in reproducing this phenomenon and therefore cannot reliably predict the density profile in the 
outer region. Therefore for r> 500 km we assume a power law p = 1.5 lO^(-R0/r)^g/cm^, which 
roughly describes the static progenitor star. 

The peculiar Fg profile in fig. |Hlis inevitably dictated by the deleptonization process |I^I15] : 
behind the shock wave which has passed in the mantle matter, the electron capture on the 
newly liberated protons is rapid, driving to low values (~ 1/4). In the outer region, where 
the density is sensibly lower, the efficiency of the capture is much lower, so that Yg essentially 
maintains the value ~ 1/2 typical of normal matter. This is important because the matter 
potential Ve of electron (anti)neutrinos flips sign, see eq. ()20p. when, in the deep region of the 
mantle Yg steeply decreases below 1/3. On the contrary, at this point the matter potentials of 
muon and tau (anti)neutrinos are marginally affected. Both are plotted in fig. |S| 

The data refer to ~0.3 sec after bounce for a typical star of ~ 11 solar masses. The subsequent 
evolution is supposed to move the wave of the profile slightly outwards, maintaining, however, 
its characteristic shape. The slight dependence on the progenitor mass, in turn, is not really 
relevant [HB] . 

This SN density profile has been computed jSOl in absence of sterile neutrinos effects. Adding 
a sterile neutrino, z/g — > Ug conversions can reduce Yg due to a non-trivial feedback mechanism on 
the MSW potential experienced by neutrinos |^ 133] . This could even create a new intermediate 
region with Y^ < 1/3, thus introducing two more level-crossing in the Ue channel. We neglect 
these possible extra MSW resonances because they do not affect the I'e rate when they are both 
adiabatic, or both fully non adiabatic. 

In summary, although the profiles that we adopt come from a specific computation and 
refer to a specific instant in time, they incorporate the peculiar features that are important for 
our purposes. A more refined treatment of this point (later times behavior of the profiles, fine 
structures connected with the passage of the shock wave...) would of course require to obtain 
first a complete simulation of the SN evolution, including the explosion. This could be needed to 
describe the signal from a future supernova, but seems unnecessarily complicated for SN1987A. 

5.2 Results 

Active/sterile mixing significantly affects SN neutrinos for 6's ~ 1. Due to MSW resonances, 
significant effects can also be present for small mixing, i.e. — ^ and 6's 7r/2. The re- 

^^Wc thank Adam Burrows for having provided us with the data, and for useful discussions. 
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cently established active/active mixings are not taken into account in older studies [531 IM] (see 
however To understand the main features it is useful to look at the pattern of possible 

level crossings, qualitatively depicted in fig. at pageUHl The three mostly active anti-neutrino 
eigenstates depend on the radius r as dictated by their measured masses and mixings (we assume 
^13 = and normal hierarchy) and by the predicted SN density profiles, and are represented in 
fig. by the three colored curves. The mostly sterile neutrino is generically represented by an 
almost-horizontal line, which is plotted in fig. EId in the specific case of small z/g/ i^s mixing and 
Am^4 ^ Aml^^. There are three possible kinds of active/sterile MSW resonances |nilE2lE31Ell: 

1. The mostly eigenstate crosses the mostly z/g eigenstate at r ~ 100 km, where Vg flips 
sign. At this point matter effects dominate over active neutrino masses, so that active 
mass eigenstates coincide with flavour eigenstates. Since Ve flips sign in a steep way this 
resonance is effective only if Aml^ > 10~^^°eV^ (different SN simulations gives values in 
this range). 

2. If the mostly sterile eigenstate is the lightest one (in our parameterization this needs 
6's>vr/4) the two eigenstates in 1. cross again at larger r. Pictorially, this second reso- 
nance is present when the sterile black line is lower than what assumed in fig. Eb- This 
MSW resonance occurs at large r where Ve is smooth, so that it is effective down to 
Am^4 > lO^^^^eV^. Again, the significant uncertainty is due to uncertainties on the SN 
density gradient. 

3. If instead the mostly sterile eigenstate is the heaviest or the next-to-heaviest state, it crosses 
one or both of two mostly u^^r eigenstates. This is the case illustrated in fig.Et". The values 
of Am24 and Amg^ determine at which r these crossings takes place, and consequently the 
flavour composition of the mostly active states at the resonance. Entering in the SN, the 
small z/g component of disappears as soon as Vg — dominates over Am^^j^.^^ The 
color of mass eigenstates in fig. illustrates these phenomena. In any case, active/sterile 
MSW resonances with the mostly Ufj^^r states do not significantly affect the z/g rate, see 
eq. (HZj). 

These considerations allow to understand fig.sEl where we plot our results for the reduction of 
the z/g rate due to sterile mixing. 

Let us start from fig. IHt that studies ui/ug mixing. Resonance 1 gives a sizable reduction 
in region 1 of fig. and resonance 2 gives a sizable reduction in region 2. Had we ignored 
solar mixing the maximal deficit would have been 100%, while in presence of solar oscillations 
the maximal effect is a ~ 80% deficit (see also [S2]). More precisely, in the interior of region 1 
one obtains = sin^ Osun^p^ because resonances 1 and 3 are fully adiabatic. In the interior of 
region 2 one obtains = sin^ ^sun'^'p^ ^ because resonance 2 is fully adiabatic and resonance 
1 irrelevant. Therefore, given the assumed initial fluxes, the z^g rate gets reduced slightly more 
strongly in region 2 than in region 1. 

In region 12 both resonances 1 and 2 are effective, and tend to compensate among each other: 
resonance 1 converts z/g into Us and resonance 2 reconverts into z/g. In region 3, resonances 3 
gives a 20% suppression of the z/g rate, that sharply terminates when Aml^ < Am'^^^. This is 

^•^When Vfj_ ~ Vr dominates over Am^^j^ the flavour composition varies from i/^ ± to i/^ and Vt- This 
active/active resonance happens at so high densities that is not relevant for our purposes. 
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Figure 9: Sterile effects in supernovae. The iso-contours correspond to a 10, 20, 30, 40, 50, 

60, 70 % deficit of the SN z/g total rate due to oscillations into sterile neutrinos. The deficit 
is measured with respect to the rate in absence of active/sterile oscillations and in presence 
of active/active oscillations (which reduce the no-oscillation rate by ~ 10%^. We shaded as 
disfavoured by SN1987A data regions with a deficit larger than 70%. While the qualitative pattern 
is robust, regions with MSW resonances can shift by one order of magnitude in Am^ using 
different SN density profiles, v^/vs mixing (not plotted) does not give significant effects. Fig. 
studies in detail the sample points here marked as A, B, C, and the regions 1, 2, 12, 3 are 
discussed in the text. 
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Figure 10: Sterile effects in supernovae. Distortion of the flux at sample points A, B, C. 

due to a strong suppression of tlie mostly t eigenstates, which due to solar oscillations would 
give a 20% contribution to the Ue rate (ignoring solar mixing, there would be no suppression 
of the Ue rate in region 3). This reduction of the z/^^t fluxes induced by resonances 3 could 
be better probed by measuring the NC rate (which gets a <40% reduction) and, if neutrinos 
cross the earth, by distortions of the energy spectrum. At smaller Am^ and around maximal 
mixing vacuum oscillations can reduce the rate by <50%: their effect persists down to 
Am^ ~ Ey/D ~ lO^^^eV^. The precise value depends on the distance D (we assumed D = 10 
kpc). 

The other mixing cases are understood in similar ways. We remind that our parametrization 
is discontinuous at 9^ = 7r/4: this is reflected in fig. I^b; where we consider z/2/z/s mixing. Res- 
onance 2 sharply terminates when Aml^ < Arri^^^^. Up to these differences, this case is quite 
similar to the previous one, because Pi and z/2 get strongly mixed by matter effects. 

On the contrary fs/z/g mixing (not shown) does not give a significant reduction of the z/g rate. 
Mixing with the flavour eigenstates behaves in a similar way. Namely, the v^^t/^s figure shows 
the reduction in region 3 due to resonances 3, and the Ve/i^s figure shows the reduction in region 
1 due to resonance 1. The additional feature at Am\^ ~ Arn^^^ is due to adiabatic conversion 
(for large sterile angles) between the mostly-sterile state and 1^2, that are almost degenerate in 
this condition. For larger Am\^, there is no Ve component in i/^ so that the crossing is totally 
non adiabatic, while for smaller Am\^ the two states are separated. As in the solar case, sterile 
effects persists at all values of Am\^, even if it is small. Unlike in the solar case, such effects are 
not sensitive to vacuum oscillations (i.e. Arn? ~ E^/D ~ 10^^* eV^) because in a SN neutrinos 
of all energies experience the adiabatic LMA resonance. Therefore in the lower row of fig.slHlwe 
only show results at Am^4 > 10~^eV^; nothing changes at smaller values. 

SN1987A data can be precisely compared with expectations doing an event by event fit [HI]; 
however the result strongly depends on the assumed average energies and total luminosities of 
the initial fluxes. The data do not permit detailed studies of the energy distribution. Even 
more, they cannot discriminate in a significant way between a larger total flux with smaller 
average energy and a smaller flux with higher average energy. In absence of a quantitative 
estimation of theoretical uncertainties, and in view of the doubtful aspects of SN1987A data 
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discussed previously, today we cannot derive precise constraints. Therefore we simply shaded as 
'disfavoured' regions where sterile effects reduce the rate by more than 70%. This is nothing 
more than a reasonable arbitrary choice. 

Future data will permit to know precisely the total rate of v^. events, and also its distribution 
in energy and time. The most important single observable could be the average Ey^ energy, 
or more generically the energy spectrum. Rather than showing iso-contour plots of (-Ep^) we 
describe their main features. In absence of sterile oscillations we expect (-Ep^) ~ 15MeV. Along 
the 'diagonal sides of the MSW triangles' 1 and 2 (e.g. around our sample point B) the average 
energy increases up to (-Ep^) ~ 18 MeV. Along the 'diagonal sides of the MSW triangle' 12 it can 
decrease down to (-Ep^) ~ H MeV. Vacuum oscillations can give the well known distortions of the 
spectrum, as exemplified in fig. HOC These effects seem larger than experimental and theoretical 
uncertainties. In all other cases sterile effects give a quasi-energy-independent suppression of the 
rate, and therefore negligibly affect {Ep^). Fig. IIOIA, gives an example. Of course, the average 
energy of positrons generated by VeV ne scatterings is higher than {EpJ because the cross 
section increases with energy. 

We focused on supernova antineutrinos. The analogous plots for neutrinos cannot obtained 
from our u plots by flipping tan 6*8 — > 1/ tan 6*8, because we are taking into account the effects of 
active/active oscillations. Effects in neutrinos are more similar to what happens to solar z/g. 

5.3 Hints and anomalies: supernovae 

To conclude, we list hints of anomalous effects possibly related to sterile effects in SN neutrinos: 

• In the past, the initial SN neutrino fluxes in different flavours were believed to follow an 
almost exact equipartition of the energy and have gaps in neutrino temperatures larger than 
in eq. (typical values were 13, 16, 23 MeV |1E])- Under these assumptions SN1987A 
data, which point to a lower average energy, disfavor solar oscillations with large mixing 
angle (sin^ 26'sun < 0.9 at 99% C.L. [HS]) while solar data established the relatively large 
mixing sin^ 29sun ~ 0.8. The tension strongly depends on the initial temperatures, and 
disappears assuming the more recent values in eq. (fT^ 

• Pulsar motion. It has been proposed that the resonant (or non resonant [0^1) conversion of 
sterile neutrinos could explain the observed large proper velocity of newly born neutrons 
stars [HZl- Indeed, in presence of very strong and axially oriented magnetic fields, which 
could be plausible in the NS environment, the MSW potential includes a contribution (at 
one loop, from neutrino scattering on the polarized medium) which depends on the relative 
angle between the neutrino momentum and the magnetic field; this leads to an asymmetric 
neutrino emission, which could be enough to account for the observed velocities. This 
mechanism needs sterile neutrinos with keV-scale masses and ^ 10"''. 

• r-process nucleosynthesis z/g Ug conversions in the mantle of the star could drive 
Ye to low values and have thus been proposed to provide a favorable environment for 
r-process nucleosynthesis. The relevant range of z/g/i^s mixing parameters spans Am^ ~ 
(1 -7- 10^) eV^ and sin^ 29s ~ 10~^ -r- fewlO"^, which should be checked against the bounds 
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from cosmological probes that we discussed in the present paper. In the same region, the 
"conversion plus re-conversion" of antineutrinos instead guarantees that the SN Vf, is not 
dramatically affected. 

5.4 Neutrinos from other astrophysical sources 

Different kinds of experiments will try to detect neutrinos emitted by extragalactic sources, such 
as active galactic nuclei. There are no firm expectations. In particular we do not know if the 
fluxes of these neutrinos will be detectably large, and eventually in which energy range. However, 
since these neutrinos are presumably mostly generated by vr decays, like atmospheric neutrinos, 
one expects similar fluxes of v and v with flavour ratio e : : r ~ 1 : 2 : (this expectation 
might be wrong). Atmospheric oscillations then convert the flavour ratio into 1:1:1, which is 
blind to other active/active oscillations. 

In presence of an extra sterile neutrino, the e:/i:r:s~l:l:l:0 flavour ratio is 
not blind to extra active/sterile oscillations. Sterile oscillations can reduce the fluxes of active 
neutrinos, and vary the relative flavour proportion. Since the total initial fluxes and energy 
spectra are unknown, and since experiments will probably be able of tagging /i and maybe r 
neutrinos well [HH]; we focus only on the observable $^/$^. 

We assume that the baseline, about 100 Mpc, is much longer than all oscillation lengths. In 
such a limit the oscillation probability reduces to multiplication of probabilities (rather than of 
quantum amplitudes): 



One can verify that sterile mixing with a flavour eigenstate ui {i = e or fi or r) mainly gives a 
depletion of $£. However, as discussed in the rest of this paper, in such a case large active/sterile 
mixing angles are already disfavoured by other data, so that it is not possible to get a sizable 
effect. 

Present data allow large active/sterile mixing with mass eigenstates i^j (z = 1 or 2 or 3) if Us 
and Ui are quasi-degenerate, with Am? <C 10~^ eV^. In view of the long base-line, neutrinos from 
cosmic sources are affected by oscillations with Am^ > 10~^^eV^. However each of the mostly 
active mass eigenstates contains roughly equal component of i/^ and Ur (unless the atmospheric 
mixing angle is non-maximal, or unless ^13 and the CP-phase signiflcantly differ from zero): 
therefore z/g/i/j oscillations do not signiflcantly affect $^/<I>t-. 

In conclusion, near-future experiments that will try to discover neutrinos from cosmic sources 
do not seem to allow promising searches of sterile neutrino oscillations. See also |69j . 

Eq. also allows to study vacuum oscillation effects in other kinds of cosmological neu- 
trinos: 

• In a near future it seems possible to detect the Ue emitted by past core-collapse su- 
pernovae. At the moment we only have order-of-magnitude predictions, but future SN 
experiments and studies might allow to better predict their total flux or spectrum. SK 
almost reached the apparently necessary sensitivity, and it seems possible to improve the 
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efficiency of tagging neutrons emitted in UeP — ^ en scatterings [70]. Relic SN Ue are af- 
fected by oscillations in the SN (as discussed in the previous section) and by oscillations 
in vacuum (down to Am^ > 10^^^ eV^). 

• In a far future, it might be possible to directly detect CMB neutrinos by coherent 
scatterings. This would allow to probe oscillations down to Am^ ~ 10~^°eV^. 

6 Sterile effects in atmospheric, reactor and beam neu- 
trinos 

In this section we discuss how SK, K2K, MACRO, Chooz, Bugey, CDHS, CCFR, Karmen, 
Nomad, Chorus and future experiments of these kinds probe sterile oscillations. KamLAND 
data have been studied in section IH together with solar data. 

6.1 Technical details 

It is convenient to compute oscillation probabilities using the neutrino and anti-neutrino density 
matrices. We convert the initial value (e.g. p{Ey) = diag ($e(-E'i/), ^/^(-Ejy), 0, 0) for atmospheric 
neutrinos) to the instantaneous mass eigenstate basis, pm = pV . In this basis evolution in each 
medium (air, mantle, core) is simply given by a diagonal matrix of phases, ^ = diag exp(— 2z(5) 
where 5i = m1^,L/ AEy. Using the matrix density in the mass eigenstate basis, we can analytically 
average 'fast' oscillations to their mean value by appropriately inserting a small imaginary part 
e in the oscillation phases 5. This is achieved by evolving the matrix density as p^{L) = 
p^(0)e^*''*^~'^''''^' where 5ij = Si — 5j and e is an arbitrary positive small number (we choose 
e = 0.01|(5jj|). This makes computations much faster than usual techniques which require lengthy 
numerical averages of the oscillation factors. In the simplest case of vacuum oscillations of 2 
neutrinos this amounts to modify the oscillation factor as sin^ S12 — > (1 — e^'"'^^^' cos25i2)/2. 

At the air/mantle and mantle/core boundaries eigenstates change in a non-adiabatic way: 
this effect is accounted by the 'level-crossing' flavour matrices P already described at pagelT^ 

Our fit of atmospheric data takes into account all most recent results: SK atmospheric 
data [7T], the K2K spectrum and total rate [72], and MACRO [72] data about trough-going 
muons. These events arise from neutrino scatterings in the rock below the detector: MACRO 
is competitive with SK because the important parameter is the surface of the detector, rather 
than its mass. We do not include data from older atmospheric experiments, which studied 
essentially the same kind of observables better measured by SK. Data are fitted by forming a 
global x^- Uncertainties are taken into account following [7J, and systematically working in 
Gaussian approximation. More naive definitions of the (^-g- fitting only the zenith-angle 
spectra of the single classes of events) do not give significantly different final results. 

Atmospheric data contain the evidence for the atmospheric anomaly. We assume it is due 
to Vj. oscillations and take into account the uncertainty on Arnl^^ and on ^atm using the 

same technique developed in section W7]\ for marginalizing over Am^^j^ and on 9^^^^. Even in the 
atmospheric case, this technique now gives an accurate analytical approximation to the usual 
active-only atmospheric fit (not shown). As previously discussed, we assume ^13 = 0. 
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Finally, we include Ve disappearance data from the Chooz [73] (14 bins) and BuGEY jTJ] 
(60 bins) reactor experiments. We include disappearance data from CDHS [70] (15 bins) 
and CCFR [77] (15 bins). NuTeV data could give additional information, if certain anomalous 
features will be understood [TS]. Although disappearance experiments give the dominant con- 
straint on sterile effects, we also include v^. data from Karmen [TU], using the likelihood 
computed by the Karmen collaboration on an event-by-event basis. ^'^ LSND results are not 
included and separately discussed in section Finally, we include Nomad and Chorus data 
on Uj. [HU], which are relevant in the case of ly-i/iys mixing at Am^ ~ 10 eV^. 

6.2 Results 

Since there are many relevant experiments, it is useful to divide them into two classes, and 
separately study their impact: 

1) Fig.s ^2 show the constraints from experiments that are not sensitive to the atmospheric 
anomaly (Chooz, Bugey, CDHS, CCFR, Karmen, Nomad, Chorus). Sterile oscilla- 
tion effects in such experiments are simply described by eq. Q, so that we only had to 
compile their results. Disappearance experiments provide the dominant constraints. 

2) Atmospheric neutrinos are a powerful probe of neutrino oscillations with large mixing 
angle and Am^ > 10~^eV^. Fig.s IT^ show how the experiments that see the atmospheric 
anomaly (SK, MACRO and K2K) constrain extra oscillations into sterile neutrinos. 

Analyzing the second class of experiments requires a non trivial work: these experiments probe 
sterile mixing in a significant but indirect way. 

The upper limit on the 'sterile fraction involved in atmospheric oscillations' conventionally 
quoted in the literature corresponds to the following 'minimal' configuration: Vs/^t mixing with 
Am34 ^ Am^jj^. In this limit, the 'sterile fraction' is related to our 6^ angle by rj^ = sin^ 9s 
and our analysis reasonably agrees with previous computations. Atmospheric neutrinos oscillate 
at a single detectable frequency. Matter effects become relevant at higher neutrino energies 
and suppress active/sterile mixing without affecting oscillations. This effect allows 

to indirectly discriminate the two channels. Sterile oscillations are mostly disfavoured by the 
zenith-angle spectra of /i-like events with TeV-scale energies. 

More general active/sterile oscillation schemes manifest in different ways. A sterile neutrino 
mixed with different flavours can give appearance signals. Active/sterile oscillations with Am^ 
comparable to Am^^j^ give oscillations at multiple frequencies, distorted by matter effects. When 
Am|4 < matter effects increase the oscillation length and give a non-standard energy and path- 
length dependence of Piv^ — > v^^)- SK can best probe such effects with its multi-GeV /i-like 
sample. However, SK cannot safely test if an oscillation dip is present. This kind of studies 
needs a dedicated future detector [HTj . 

The SK collaboration probes sterile effects in another more direct way: selecting a sample 
of NC-enriched events [TJ. In absence of a precise description of the cuts performed to obtain 
the NC-enrichment, we have not included this sample (as also done in previous reanalyses of SK 
data). The 'minimal' active/sterile configuration is more strongly constrained indirectly by the 

^''Wc thank K. Eitcl and M. Steidl for giving us the tabic of the Karmen hkehhood, and B. Louis and G. Mills 
for the LSND hkclihood. 
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Figure 11: Sterile mixing: effects in short base-line experiments Chooz, Bugey, 
CDHS, CCFR, Karmen, Nomad, Chorus. Shaded regions: excluded at 90,99% C.L. The 
blue dashed lines estimate the region that seems explorable by a future short-baseline reactor 
experiment. The plot is symmetric under tanOs ^ 1/ tan6's so that we only show tanOs < 1. 



zenith-angle spectra. We estimate that this remains true in most of the parameter space, and 
consequently we do not include in our final results our approximate reanalysis of NC-enriched 
data. 

As in the solar case, we looked if present data contain some evidence for sterile effects which 
correct in a minor way many observables, by searching for local minima of the global x^- No 
statistically significant hint is found: since subleading sterile effects do not improve the global fit 
in a significant way (at most by Ax^ ~ 4) our plots only show excluded regions. The excluded 
region in fig.sll2li. e (which correspond to u^/ Ug and to Ur/ Vs mixing) extends down to Am\^ = 
because even in this limit there are sterile oscillations at the atmospheric and solar frequencies. 

It is useful to compare the sensitivity of the two classes of experiments, 1) and 2). Since 
there are no MSW resonances, all these experiments are sensitive only to relatively large sterile 
mixing, ^s>0.1. Sterile mixing with i/g (and with the ui and z/2 mass eigenstates that contain 
a sizable z/e fraction) is better probed by reactor experiments, although e-like events at SK 
extend the sensitivity down to smaller values of Am^. On the contrary atmospheric experiments 
give more stringent tests of Ug/ Vr mixing and of z/g/ mixing. Within standard cosmology, the 
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Figure 12: Sterile mixing: effects in atmospheric neutrinos (SK, MACRO) and K2K. 

No statistically significant evidence is found. Shaded regions: excluded at 90, 99% C.L. Coloured 
lines are iso-curves of a few promising signals. Blue dashed lines: 5% and 1% reduction of the 
NC rate at MINOS. Green dot-dashed lines: Pe^ = 0.01. Red dotted lines: |AP^t-| = 0.01 at 
CNGS. 

sterile effects detectable by the experiments discussed in this section are already disfavoured by 
measurements of the primordial ^He abundancy, and can be fully tested by future CMB or BBN 
data. 

A detailed analysis of capabilities of future beam or reactor neutrinos as probes of sterile neu- 
trinos seems not necessary. In fact, there are many proposals motivated by other considerations, 
and in each case it is easy to compute sterile effects. We only make a few general comments. 

The blue dashed line in fig.s ITTl shows what can achieved by a future high-precision short- 
baseline reactor experiment able of detecting a 2% deficit in the Ve flux [H2]- Sterile oscillations 
give a Ve deficit, which might be energy-dependent if the sterile oscillation length at Ey ~ 
fewMeV is comparable to the base-line L. Both are unknown; we assumed L ~ 2 km. Most of 
the explorable region at small Am? is already excluded by solar and atmospheric experiments. 
The region with large Am^ ~ eV^ is more difficult because even a near detector only sees 
averaged oscillations; one has to rely on theoretical predictions for the total flux of reactor 

i/e IBS!. 
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Sterile oscillations with large Am^ ~ eV^ have a wave-length comparable to the earth radius 
at energies E^, ~ TeV: experiments such as AUGER and IceCube can study atmospheric 
neutrinos in this energy range, and could see an anomalous zenith-angle dependence if is large 
enough [M] , 

Future atmospheric experiment such as Monolith jH^ could test if the first oscillation dip is 
present as predicted by z/^ — > oscillations. Sterile effects with Am^ ~ Am'^^^ can significantly 
distort the expected oscillation pattern. 

Distortions not related to earth matter effects can also be searched by the planned z/^-beam 
experiments. Furthermore, the MINOS detector can statistically distinguish NC/CC events 
that are tagged as short/ long tracks, and it should be possible to measure the NC rate with 5% 
precision jHS]- The blue dashed lines in fig.s El correspond to a 5% and 1% anomaly in the NC 
total rate in a z/^ beam experiment with L = 730 km and MiNOS-like parent energy spectrum. 

A sterile neutrino can also manifest as u^- appearance: considering a CNGS-like beam [Sni in 
fig.s^lwe show iso-curves corresponding to an average z/^ — conversion probability ±0.01 
different from atmospheric oscillations only (which give a conversion probability of about 0.02, 
depending on the precise value of Aml^^ and on the average energy of the beam). 

Constraints from reactor experiments make more difficult, but not impossible, to have de- 
tectable i>e appearance in a z/^ beam (or t'^ appearance in a Ue beam) as signals of active/sterile 
oscillations. In experiments sensitive to atmospheric oscillations a non zero 6is gives such ef- 
fects. Assuming ^13 = 0, in fig. El we show iso-curves corresponding to a z/^ — z/g conversion 
probability of 0.01. 

6.3 Hints and anomalies 

To conclude, we list some possibly anomalous features of present data: 

• The total rate of e-like events at SK is higher than the expected central value. The n/e 
ratio seems to be too low (see e.g. |HZ|)- Constraints from solar, atmospheric and short- 
baseline data show that sterile neutrinos can reduce the /i rate (e.g. in the case of t'^/z/g 
mixing with Amf^ ~ 10~^eV^) but cannot increase the e rate. 

• The total rate of trough-going muons at SK is higher than the expected central value. 
Predictions and direct measurements of cosmic ray primaries at the relevant energies are 
difficult IHHl. 

• The LSND result jl], discussed below. 

6.4 LSND 

The LSND experiment |1| finds an evidence for ^ Ve-, that ranges between 3 to 7cr depending 
on how data are analyzed. It is difficult to conceive new physics that can explain the LSND 
result compatibly with all other constraints. Among reasonably conservative interpretations, 
a sterile neutrino with eV-scale mass seems to be the less disfavoured possibility Ej- 
According to this interpretation the LSND anomaly arises as t'/i — »• z/g — > z/g so that the effective 
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Figure 13: The LSND anomaly interpreted as oscillations of 3+1 neutrinos. Shaded 
region: suggested at 99% C.L. by LSND. Black dotted line: 99% C.L. global constraint from 
other neutrino experiments (mainly Karmen, Bugey, SK, CDHS). Continuos red line: N^, = 3.8 
thermalized neutrinos. Dot-dashed orange line: Vt^h'^ = 0.01. 

^LSND ^fi/^e mixing angle is predicted to be ^lsnd ~ ^es ■ 0^J.s■ This formula is valid only for small 
mixing angles, and eq. (j21) gives the general expression. This prediction gives rise to 3 problems: 

1) i/g and disappearance experiments imply that 6es are 6^s are somewhat smaller than 
what suggested by LSND [HniE2lEl; 

2) standard BBN predicts that the sterile neutrino thermalizes, so that primordial abundances 
should have values corresponding to A^^, = 4 PU] : 

3) according to standard cosmology, the sterile neutrino gives a contribution to the neutrino 
density fli, somewhat larger than what suggested by global fit of CMB and LSS data (see 
e.g. ED). 

Concerning points 2) and 3), there is not yet general consensus that the sterile neutrino ther- 
malizes, maybe because this LSND issue has never been analyzed by authors that performed, 
at the same time, a precise study of neutrino data and of cosmology with mixed neutrinos. Es- 
timates indicate that the region favoured by LSND lies well inside the region where the sterile 
neutrino is thermalized jHOl- This is confirmed by our analysis, shown in fig. E| In the rele- 
vant region the constraint on Q^y is well approximated by the horizontal line corresponding to 
fl^h"^ = m4/93.5 eV, as assumed in previous analyses jnUE]- The accurately computed constraint 
starts to be weaker only at much smaller values of the effective 6'LgND mixing angle. 

In fig. ^1 the BBN constraint has been minimized (when allowed by neutrino data) setting 

1 /2 

0es ~ ^/is ~ ^LSND- We scc that short-baseline experiments sensible to a P{i^^ I'e) about 2 
orders of magnitude smaller than the value suggested by LSND are needed to probe regions 
compatible with standard BBN. 
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7 Summary 



A few years ago active/sterile oscillations were studied as an alternative to active/active oscil- 
lations. For example, it was shown that neither the solar nor the atmospheric anomaly can 
be produced by oscillations into sterile neutrinos compatibly with standard BBN. These active 
vs sterile issues have now been firmly solved by experiments, and the new relevant questions 
become: 

How large can be the subdominant sterile component possibly present in solar or 
atmospheric oscillations? How can we discover new anomalies due to sterile effects? 

In order to address these issues we systematically compared oscillation effects generated by one 
sterile neutrino (including the effects of the now established solar and atmospheric oscillations) 
with present experiments and studied capabilities of future probes, extending previous analyses 
in several ways. Almost all previous analyses that studied these new questions considered only 
the peculiar sterile oscillation pattern that gives the simplest physics: the sterile neutrino mixes 
with a specific flavour (z/^^,- in solar analyses, Vr in atmospheric ones) and has a large mass, 
Am^ ^ Am^^j^gy^. Dropping each one of these assumptions gives quite different physics. We 
explored the whole parameter space (three sterile mixing parameters and one sterile mass) 
producing precise results for six representative slices, that span the spectrum of the various 
possibilities. We considered one extra sterile neutrino with arbitrary mass m4, and allowed it to 
mix with 

Ve oi or or Vi or V2 01 

where 1^1,2,3 are the mass eigenstates in absence of sterile mixing. The spectrum of active neutri- 
nos is not yet fully known: rather than studying all possible cases we focussed on what we believe 
is most plausible case: we assumed normal hierarchy of active neutrinos (i.e. mi <^ 7712 ~ rris/G) 
and that ^13 is small enough that we can neglect its effects. If experiments will contradict 
these assumptions, it will be easy to update our results; most of them do not depend on these 
assumptions. 

We considered the most promising ways to probe the existence of eV-scale sterile neutrinos. 
Most probes are based on a careful study of natural sources of neutrinos (the universe, the 
sun, supernovae, cosmic rays,...) which have their own peculiar capabilities and limitations. We 
studied how one extra sterile neutrino affects BBN (helium-4, deuterium), CMB, LSS and solar, 
atmospheric, reactor, beam experiments. The sensitivity of some of these probes is enhanced 
by MSW resonances [221 • In cosmology, active u and u encounter a MSW resonance with 
sterile neutrinos lighter than active ones. Roughly the same happens to supernova Ue (that also 
experience less important MSW resonances in the opposite situation). On the contrary, solar Uf, 
encounter a MSW resonance with sterile neutrinos heavier than active ones. 

Fig.s El combines present constraints. Each probe is described in greater detail in its specific 
series of figures: fig.sElfor cosmology, fig.slHlfor solar experiments, fig.slHlfor supernovas, fig.sITT] 
for short base-line experiments, fig.s ^1 for atmospheric experiments. These figures also show 
the capabilities of some future experiments, that we now try to summarize in words. 

• Compatibility with standard BBN constrains sterile oscillations occurred at temperatures 
T> O.lMeV. It is very important to improve measurements of the helium-4 primordial 
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Figure 14: Summary of sterile neutrino effects. The shaded region is excluded at 99% 
C.L. (2 dof) by solar or atmospheric or reactor or short base-line experiments. We shaded as 
excluded also regions where sterile neutrinos suppress the SN1987A Ue rate by more than 70%. 
This rate is suppressed by more than 20% inside the dashed blue line, that can be explored at the 
next SN explosion if it will be possible to understand the collapse well enough. Within standard 
cosmology, the region above the red continuous line is disfavoured (maybe already excluded) by 
BBN and LSS. Plots in the various sections show how much each probe can be improved by future 
experiments. 
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abundancy (that we parameterize in terms of an effective number of neutrinos Nj^*^, see 
eq. (jU) until A'^^,^'^ = 4 will be safely tested. This requires overcoming 'systematic' uncer- 
tainties. The helium-4 abundancy is sensitive to two different sterile effects: increase of the 
total neutrino density, and depletion of electron- neutrinos. The second effect makes N^^^'^ 
sensitive to sterile oscillations down to Am^ ~ 10~^eV^, while the first effect becomes 
negligible at Am^ < 10~^eV^. 

However, if BBN were non-standard, a modified density of electron neutrinos could com- 
pensate the sterile corrections to A'^^^'': for example helium-4 constraints on sterile oscil- 
lations can be evaded by allowing a large neutrino asymmetry. 

For all these reasons it is important to measure a second BBN effect. The deuterium 
primordial abundancy is affected by milder systematic problems: in the future it might be 
possible to improve its measurement obtaining an uncertainty on the effective parameter 
A'^^ (precisely defined in eq. dH))) significantly below 1, possibly making deuterium the 
most significant BBN probe. We have computed the ranges of active/sterile oscillation 
parameters that significantly affect N^: it is less sensitive than helium-4 to depletion 
and therefore to values of Am^ below 10~^eV^ (fig- El)- 

Future studies of Cosmic Microwave Background acoustic oscillations should allow to 
precisely measure the total neutrino density N^^^ at recombination (T ~ eV) with ibO.2 
(Planck) or maybe ±0.05 (CMBpol) error [T7j. Neutrinos affect CMB in various ways; 
neutrino free-streaming offers a clean signature that allows to count neutrinos. Sterile 
neutrinos affect N^^^ only if Am2>10-^eV^; in such a case A^^^^ ^ A^^^° ^ N^. 
Therefore CMB will cover only a part of the region that BBN could probe with fully 
reliable measurements of the helium-4 or deuterium abundances. 

CMB also allows to probe eV-scale u masses. Smaller u masses can be probed by mea- 
suring how much galaxies are clustered, because neutrinos become non relativistic when 
the observable universe had a size comparable to present cluster of galaxies. Relativistic 
neutrinos freely move and tend to reduce the amount of clustering. Recently, Large Scale 
Structure data (together with precise CMB measurements) gave a bound on the present 
energy density in neutrinos < 0.0076 at 95 % C.L. [SI HH], dominated by neutrino 
masses (rather than neutrino energy). With only active neutrinos this implies a significant 
bound on their mass, rrii, < 0.23 eV [51IT^. 

Sterile neutrinos contribute to fl,^ (a small non-thermal population of relatively heavy 
sterile neutrinos also modifies the way fl^, manifests). For small active/sterile mixing LSS 
constrains sterile oscillations more strongly than BBN (fig. E}. 

Solar Ue experiments have explored sterile oscillations not testable by BBN, thanks to 
two different effects. (1) MSW resonances make solar z/g sensitive to small active/sterile 
mixme; and Am2 > 10"^ eVl (2) With large mixing, solar are sensitive down to Am^ ~ 
10~^^eV^. Future experiments will explore new aspects of the solar neutrino anomaly, 
allowing to measure in a redundant way the active oscillation parameters or to discover 
a new anomaly. We emphasize one qualitative point. Due to LMA oscillations, neutrinos 
exit from the sun as almost pure V2 at energies > few MeV. Neutrinos with these 
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energies have been precisely studied by SK and SNO, but are almost unaffected by sterile 
oscillations if they involve mostly z/i. This could happen either when z/g mixes with i>i, 
or when Us experiences a level-crossing with ui. Therefore there is a whole class of sterile 
effects which manifest only at < few MeV — an energy range explored so far only by 
Gallium experiments. Future precise measurement of solar at sub-MeV energies will 
allow to significantly extend searches for active/sterile effects. Part of these extended 
region can be soon tested by Borexino, where a sterile neutrino can manifest as day/night 
variations, or as seasonal variations, or even by reducing the total rate. 

• Supernova neutrinos will be good probes of sterile oscillations because have a different 
pattern of MSW resonances and a longer base-line than solar u^. Consequently supernova 
Ue are more sensitive than solar in two main cases: (a) small Am^>10~^^eV^ with 
large Og] (b) Ug lighter than ui with small mixing. Oscillations into one sterile neutrino can 
reduce the z/g rate by up to 80% (see fig.sEl) ^^d, in a more restricted range of oscillation 
parameters, vary the average z/g energy by 30%. SN1987A data agreed with expectations. 
Future SN experiments can perform quantitative test, but it is not clear how to deal with 
theoretical uncertainties. We also discussed other less promising astrophysical probes. 

• Atmospheric experiments (SK, MACRO, K2K) indirectly exclude active/sterile oscilla- 
tions with Am? > 10~^"^^ eV^ and tan^ ^ 0.1 ^ 0.2. Up to minor differences, this applies 
to all flavours (fig. fT2|l . Terrestrial experiments that mainly probed disappearance of 
z/g and z/^ (Chooz, CDHS,...) exclude active/sterile mixings with these flavours with 
tan^^s^O.03 and Am^>10~'^eV^ (fig. ITTjl. Therefore future short-baseline experiments 
can search for sterile effects with smaller 6s- Possible signals are disappearance in reactor 
experiments, a deficit of NC events or u^- appearance in beam experiments. Within stan- 
dard cosmology these effects can be probed by CMB and BEN, which already disfavour 
them. 

We listed present anomalies that can be interpreted as due to sterile neutrinos. None looks 
particularly significant, with the possible exception of the LSND anomaly, that can be due 
to oscillations of '3 + 1' neutrinos. Our precise study confirms that the extra sterile neutrino 
suggested by LSND thermalizes almost completely before BEN (in agreement with estimations 
of most previous analyses). In fig. we plot, as function of the effective LSND oscillation 
parameters ^lsnd and AmLg^D, the constraints on A^^^*^ ~ A^^ ~ A^^^^ (mainly from helium-4 
data) and Q,^ (mainly from LSS data): in the relevant region our precise plot negligibly differs 
from previous estimates [Hllini- The 3+1 interpretation of LSND is not compatible with standard 
EEN, and gives a which is only marginally compatible with standard cosmology. 
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